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Abstract 



In this paper we study the splitting of separatrices phenomenon which arises when one considers a 
Hamiltonian System of one degree of freedom with a fast periodic or quasiperiodic and meromorphic 
in the state variables perturbation. The obtained results are different from the previous ones in 
the literature, which mainly assume algebraic or trigonometric polynomial dependence on the state 
^ \ variables. As a model, we consider the pendulum equation with several meromorphic perturbations 

^^ ■ and we show the sensitivity of the size of the splitting on the width of the analyticity strip of the 

perturbation with respect to the state variables. We show that the size of the splitting is exponentially 
small if the strip of analyticity is wide enough. Furthermore, we see that the splitting grows as the 
width of the analyticity strip shrinks, even becoming non-exponentially small for very narrow strips. 

C ' Our results prevent from using polynomial truncations of the meromorphic perturbation to compute 

^— ^ , the size of the splitting of separatrices. 



1 Introduction 



. , Exponentially small splitting of separatrices appears in analytic dynamical systems with different time 

j^ ' scales. A paradigmatic example are analytic Hamiltonian systems of one degree of freedom with a fast 

non-autonomous periodic or quasiperiodic perturbation. Namely, systems of the form 

H {x,y,-\ = Ho{x, y) + fj,e''Hi ix,y,- 

= y- + v{x)+fie'^H,(x,y,^-Y 

where e > is a small parameter, /i G M, 77 > 0, and Hq has a hyperbolic critical point whose invariant 
manifolds coincide along a separatrix. 
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This phenomenon was first pointed out by Poincare |Poi99] but it was not until the last decades when 
this pr oblem started to be studied r igorousl y (see for instance fHMS88, SMH9f. DS92. Fon93l [CG9l 
[Gil9llFto95l ltJau95[ inUJ597l 105971 [Gil97bl [I7i97l IG!GM99I iGilOO. .SauOl. ,DGS04. ,BF04. ,BF05l iBilOel 
IOli06[ [GOSlOl iBFGSll] ). Nevertheless, the results which obtain asymptotic formulas for the splitting 
only deal with Hamiltonian systems whose perturbation is an algebraic or trigonometric polynomial with 
respect to the state variables x and y. In all these cases the splitting of separatrices is exponentially 
small with respect to the parameter e. Moreover, the imaginary part of the complex singularity of the 
time-parameterization of the unperturbed separatrix closest to the real axis plays a significant role. 

All the previous works dealing with the periodic case, show that under certain non-degeneracy con- 
ditions, the distance between the invariant manifolds is of order 

d-Me^e't, (2) 

where a is the imaginary part of the complex singularity of the time-parameterization of the unperturbed 
separatrix closest to the real axis and q G M. Moreover, for rj > rj* , where 77* depends on the properties of 
both Hq and Hi, the splitting is well predicted by the Poincare- Arnol'd-Melnikov method (see |Mel63j . 
and JGH83J for a more modern exposition of this method). This case is usually called regular case. In the 
singular case tj = rf the splitting is exponentially small as ^ but the first order does not coincide with 
the Melnikov prediction (see [BFGSll] and references therein). In the quasiperiodic case, under certain 
hypotheses, one can also show that the Melnikov method predicts correctly the splitting and the size of 
both the Melnikov function and the splitting depends strongly on a |DG JS97"1 ISauOli IDGS04J . However, 
this case is much less understood and there are very few results. 

Nevertheless, many of the models known, for instance in celestial mechanics, are not algebraic or 
trigonometric polynomials in the state variables but involve functions with a finite strip of analyticity 
(see, for instance, jLSSOi IXiaMl IMP941 IFGKRllp . As far as the authors know, the only result deahng 
with the exponentially small splitting of separatrices in the periodic case for non-entire perturbations is 
|Gel97a| . However, the author considers models with a strip of analyticity very big with respect to e so 
that he can deal with them as if they were polynomial. In the quasiperiodic case, as far as the authors 
know, there are not previous results. 

The goal of this paper is to study how the splitting of separatrices behavior depends on the width of 
the analyticity strip when one considers a meromorphic perturbation. Essentially, we see that the size of 
the splitting depends strongly on this width and that, in general, the singularity of the separatrix does 
not play any role in this size. We consider also the case when the strip tends to infinity as e — )■ and we 
see how the size of the splitting tends to the size known for the entire cases. In the other limiting case, 
namely when the strip of analyticity shrinks to the real line as £ — > 0, we see that even if the perturbation 
is still analytic, the splitting becomes algebraic in e both in the periodic and the quasiperiodic case. 

We focus our study in particular examples, which allow us to analyze in great detail the behavior of 
the splitting. Nevertheless, we expect the same to happen for fairly general systems. 

Wc work with time periodic and quasiperiodic perturbations of the classical pendulum. More con- 
cretely, we consider the following model, 

* = ^^""+^'"(rT^^KO' ^^^ 

where /(t) is an analytic function which depends either periodically or quasiperiodically on r. 
The associated system 

^ x = y 

. . sinx Jt\ (4) 

y = sm.T + l-ie'—- -. rj/ - 

(1 + asmx)"^ \e^ 

is Hamiltonian with Hamiltonian function 



H[x,y,-\ ^^+cosx-l + fis^^{x)f{-], (5) 



where ^(x) is defined by il^' {x) ~ — sinx'/(l + asinx)^ and V'(O) = 0. Here a e [0, 1) is a parameter 
which changes the width of the analyticity strip of ^j, which is given by 



, , / 1 + Vl - a^ \ 
|Inix| < hi . (6) 

When a — 0, the system is entire in a; and y and has been previously studied for particular choices of / 
in |Tre97[ IDGJS971 IOSS03[ [QhOG] . whereas when a == 1, V' is not defined in a; == 37r/2. In this paper, we 
consider any a € (0,1) either independent of or dependent on e. 

In the periodic case, as experts know, the only important property to obtain the asymptotic formula 
for the splitting is that one has to require that the first harmonics of / are different from zero. Thus, we 
choose 

/(r) =sinT. 

Dealing with any other function with non-zero first Fourier coefficients is analogous. In this setting, one 
can rephrase system (|4]) as a Hamiltonian system of two degrees of freedom considering r — t/e as a new 
angle and / its conjugate action, which gives the Hamiltonian 

K (x, y, T, I) =- + H {x, y, r) 

= — I h cos a; — 1 + iie^'ip{x) sin r. 

e 2 

In the quasiperiodic case, we consider the same model ([3]) with /(r) = i^(r, 7r). where 

^/5 + l ,„. 

7=^— (8) 

is the golden mean number and i^ : T^ — > R. Note that if one takes a = 0, one recovers the model 
considered in JDGJS97] . On the function F we assume the same hypotheses that are assumed in that 
article. Namely, if one considers its Fourier expansion in the angles 9 = {9i,62), 

F(0i,02)- 5^i^[''le*^ (9) 

we assume that there exist constants ri , r2 > such that 

^[fe]g'-l|fcl|+'-2|fc2| 



sup 

fe=(fci,fc2)ez 



< oo. (10) 



Furthermore, we assume that there exist a and fco such that 

i^I'^1 > ae-'^il'^il-''=l'=^l (11) 

for all |fci|/|fc2| which are continuous fraction convergents of 7 and \k2\ > kg. An example of function 
satisfying these hypotheses is 

m,^^) cos^icos^^ 



(coshri — cos 6*1) (cosh r2 — cos6'2) 



Introducing the angle coordinates {01,62) and their conjugate actions (/i,/2), system @ can be seen 
as a 3 degrees of freedom Hamiltonian System with Hamiltonian 

Kix,y,0,l) = — + H{x,y,0) 



e 



(12) 



'^^^— + ^ + cosx - 1 + ns''ij{x)Fiei,e2) 

e 2 



where uj — (1,7) is the frequency vector. 

We have chosen these particular models for several reasons. First, the hyperbolic critical point (0, 0) of 
the unperturbed pendulum persists when the perturbation is added. This fact is not crucial but simplifies 
the computations. Second, with the chosen function -0, the size of the Melnikov function depends on the 
strip of analyticity of the perturbations, as is expected to happen for general systems. In Remark 12.51 in 
Section [5J we consider the non-generic model 

sinx . t 

X = sm X + lie ' — sm - , 

(1 — acosa;)^ e 

which has the same strip of analyticity ([B]). However, due to certain cancellations, the size of Melnikov 
function does not depend on this strip. 

In the quasiperiodic case, we have chosen a very specific function F. On one hand, we have chosen the 
frequency vector w = (1,7), where 7 is the golden mean ([5]). The size of the splitting strongly depends 
on the diophantine properties of the chosen frequency and, in fact, its rigorous study has been only 
done, as far as the authors know, for quadratic frequencies (see |Sau01[ ILMS031 |DG03| ). On the other 
hand, the chosen function F has finite strip of analyticity in the angles (0i,02)- This is the only kind of 
systems for which it is known that the Melnikov function predicts correctly the size of the splitting (see 
|Sim94[rSV01j ). In fact, in the quasiperiodic case, the width of the strip of analyticity of F also plays a 
crucial role in the size of the splitting. 

Finally, as we have already explained, this particular choice of the perturbation makes everything easily 
computable. This allows us to obtain explicit formulas for the first order of the splitting of separatrices, 
using the Melnikov function, and see how it depends on the width of the analyticity strip of ip. Then, we 
can compare our results for a small with the existing previous ones for a = 0. 

As we have already said, when a ~ I system @ is not defined at a; = 37r/2. Therefore, it has no 
sense to study the splitting problem for a too close to 1. Indeed, the perturbation is small in the real 
line provided 

Nevertheless, as usually happen in the exponentially small splitting problems (see |GOS10p . we will see 
that the splitting problem has sense under the slightly weaker hypothesis 

<1. 



(1 - a)3/2 

When /i = 0, the system is the classical pendulum. It has a hyperbolic critical point at (0, 0) whose 
invariant manifolds coincide along two separatrices (see Figure [T]). We focus our attention on the positive 
one, which can be parameterized as 

2 

xo(u) = 4arctan(e") , yoiu) — xq(u) ~ — , (13) 

cosh u 

whose singularities are at m = iTr/2 + ikn, fc G Z. 

The goal of this paper is to measure how this separatrix splits when one takes ^ > 0, paying special 
attention on how this splitting depends on the relative size between a and e. The previous results in the 
periodic case |Tre971 IOSS031 IOh06| consider a = 0. They obtain an asymptotic formula for the distance 
between the invariant manifolds, which is of the form 



fj.e'T'^e 2£. 



In this paper, we see that in the periodic case 



• If a < e^, that is when the strip of analyticity is of order 21n(l/e) or bigger, the distance between 
the invariant manifolds coincides with the results obtained in the case a = 0, when the perturbation 
is a trigonometric polynomial. 




Figure 1: Phase portrait of the pendulum. In it, one can see the two separatrices. We focus our study 
on the upper one. 

• If e'^ <C a and 1 — a ^ e^, that is when the strip of analyticity is between 21n(l/e) and e, the 
distance is exponentially small but the first order is given by a completely different formula 

c 

where g e R and < c < tt/2 are constants which depend on a. Therefore, the splitting is bigger 
than it was in the polynomial case. Moreover, we see that it monotonously increases with a. 

• Finally, if < 1 — a > e^, that is when the strip of analyticity is narrower than e, the distance 
between the invariant manifolds is non-exponentially small. 



where g G M is a constant which depend on a. 

In particular, the second statement shows that even if a is a small parameter, our study prevents of 
using the classical approach in perturbation theory. It consists in expanding the perturbation term in 
powers of a and studying the splitting of the first order, which has a polynomial perturbation. We will 
see that even for a ~ £, this approach leads to a wrong result for the splitting. Therefore, we show that 
the study of the splitting of separatrices for meromorphic perturbations cannot be reduced to the study 
of simplified polynomial models. In the quasiperiodic case, an analogous phenomenon happens. 

The structure of the paper goes as follows. First in Section [5] we give the main results considering 
the periodic case. In Proposition 12.11 we show the behavior of the Melnikov function with respect to 
a. In Corollary 12.61 we analyze the range of the parameter a for which the Melnikov function is not 
exponentially small. Then, in Theorems 12.81 and 12. 101 we show for which range of parameters e, a and rj, 
the Melnikov function predicts correctly the splitting. 

In Section [3] we consider the quasiperiodic case. First in Proposition 13.11 we study the size of the 
Melnikov function and in Theorem 13.51 we prove its validity. As in the periodic case, we also show that 
if 1 — a > e^ both the Melnikov function and the splitting of separatrices are not exponentially small. 
These results are given in Corollarv 13.31 and Theorem [3?7] respectively. 

Section H] is devoted to prove Theorems 12.81 and 12. 101 and Section [5] is devoted to prove Theorems 13.51 
andO 

Finally, in Appendix [^ we give some heuristic ideas about how to deal with the so-called singular 
periodic case, namely, when the parameter rj reaches a certain limiting value and therefore the Melnikov 
function docs not predict correctly the size of the splitting. 



1.1 Heuristics on the relation between regularity and Arnol'd diffusion 

Even if in this paper we study the so cahed isochronous case, where the frequency of the perturbation is 
fixed, the same kind of study would apply to an anisochronous case (see |Sau01| ). In the anisochronous 
case one can encounter both situations, the case of rationally dependent frequencies, which leads to 
a periodic in time perturbation, and the case of rationally independent frequencies, which leads to a 
quasiperiodic perturbation. 

Anisochronous systems have attracted a lot of attention because they are good models to study the 
phenomenon called Arnol'd diffusion. The name comes from the fact that was V. Arnol'd who produced 
in 1964 |Arn64| the first example showing a possible mechanism that leads to global instabilities in nearly 
integrable Hamiltonian Systems. 

Arnol'd proved the presence of instabilities in the following particular model 

P P 

H{Ii,l2,'fi,'P2,t) = Y + Y +e{cosipi - 1) + fie{cosipi - l)(sinv32 +sini), 

showing the existence of orbits whose action I2 changes drastically. 

Nevertheless, it is expected that instabilities exist in fairly general nearly integrable Hamiltonian 



Systems. Chierchia and Gallavotti in |CG94| proposed the study of the following generalization of the 
Arnol'd model 



r2 



P P 

H{Ii,l2,fi,'^2,t) = Y + Y +^(cos<pi - 1) + i-Leh{ipi,ip2,t]e). (14) 

In this setting, they coined the terminology a priori stable versus a priori unstable. A priori stable refers 
to consider fi = e^ with rj>Q and a priori unstable refers to e = 1 and /i small. In the a priori stable case 
the unperturbed system, e = is completely integrable in the sense that it is written in global action- 
angle variables. In the a priori unstable case the unperturbed system, ju = 0, even if it is integrable in 
the sense that it has conserved quantities, presents some hyperbolicity, namely has partially hyperbolic 
tori with homoclinic trajectories. 

There have been some recent works proving the existence of instabilities for a priori unstable systems 
using the ideas proposed by V. Arnol'd (see |DdlLS06[ IDH09J and see }CY04[ ITre04[ IBer08| for proofs 
using other methods). One of the main steps in the proof is to see that the stable and unstable manifolds 
of the partially hyperbolic tori, which coincide when /i = 0, split producing chains of hcteroclinic orbits. 
To detect the splitting of these invariant manifolds becomes an essential step in these geometric methods. 
The main reason that makes very difficult to detect this splitting when ^ and e are small is that one 
expects this splitting to be exponentially small in e and, thus, very difficult to study. This was the origin 
of the distinction between a priori stable and unstable systems, Arnol'd overcame this difficulty taking 
the parameter jjl exponentially small in e. It is an open problem to see whether Arnol'd mechanism works 
in the a priori stable setting, namely taking /i = e'' with 77 > 0. 

Nonetheless, the regularity of the system plays a crucial role in the classification between a priori 
stable and unstable systems. The first observation, that is commonly accepted, is that in the C case, 
even in the a priori stable setting, the splitting is not exponentially small. Nevertheless, as far as the 
authors know, the only proof of this fact is the paper |DJSG99] where it is seen that for quasiperiodic C" 
perturbations the splitting is polynomially small in e. Therefore the distinction between a priori stable 
and unstable systems regarding the splitting problem only has sense for analytic perturbations. Moreover, 
the results stated in Theorems 12.101 and 13.71 show that for analytic perturbations with narrow strip of 
analyticity, the size of the splitting is not exponentially small. In particular, we see that the splitting is 
polynomial in e in both the periodic and the quasiperiodic case. Therefore, one would expect that for 
anisochronous systems (|14p with a narrow strip of analyticity the splitting at the resonances and in the 
nonresonant zones are of the same order and polynomial with respect to e. 

On the other hand, it is a commonly accepted fact that the bigger the size of the splitting the faster the 
diffusion. In fact, Nekhroshev type lower bounds of the diffusion time (or upper bounds of the stability 
time) increase with the regularity of the system |Nek77[ IMS021 IBoulOj and this agrees with the fact that 
the splitting decreases with the regularity. For analytic Hamiltonians, these bounds are bigger for bigger 
strips of analyticity |P6s93[ IDG96| . The results in this paper show that this is also consistent with the 
fact that the splitting decreases when the strip of analyticity increases. 
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2 The Melnikov function and its validity in the periodic case 

As we want to deal with a time periodic perturbation of the pendulum equation, the dynamics is better 
understood in the three-dimensional extended phase space {x,y,T) E T x M x T. In this space, A = 
{(0, 0, r);T G T} is a hyperbolic periodic orbit and, for ^ = 0, its 2-dimcnsional stable and unstable 
invariant manifolds coincide along the homoclinic manifold 

W«(A) = W^(A) = {{x, y, r) : i/o(x, y) - 0} = {{x, y, r) = (xo(w), yo(w), r); (u, r) e M x T} , 

where {xo{u), yoiu)) is the parameterization of the separatrix given in (|13p . 

Our goal is to study how these manifolds split when ^ ^ 0. To this end we need to introduce some 
notion of distance between them. As the manifolds are graphs for /i = 0, the same happens for /is'' > 
small enough in suitable domains. More concretely, in Section 2] we will see that one can parameterize 
the perturbed stable and unstable manifolds as 




Here, (w, r) £ (— oo,C/) x T, for certain [/ > 0, for the unstable manifold and {u,t) E (— [/, oo) x T for 
the stable one. Taking into account that the invariant manifolds are Lagrangian, in Section 2] we use the 
Hamilton- Jacobi equation to see that the functions y"'* can be given as 

yo{u) 

for certain generating functions T'"'*'. 

Therefore, a natural way to measure the difference between the manifolds is to compute 

D(u,T) = a„T^(M,T)-a„T"(u,T). 

If one considers a perturbative approach taking fi as small parameter, one can easily see that the first 
order in n of the function D{u, r) is given by the Melnikov function 

M{u,t)^ {Ho,Hi}(^xoiu + s),yo{u + s),T + ^^ds. (15) 

In other words, one has that 

D{u, t) = ^le"Miu, t) + (/i^e^") • (16) 

To see that the manifolds split, we can choose a transversal section to the unperturbed separatrix to 
measure their distance. The simplest one is a; = tt, which corresponds to compute D{0,t) (see (J13p ). 
Then, the zeros of D{0,t) correspond to homoclinic orbits of the perturbed system and the distance 
between the invariant manifolds is given by 

dir) = y^(0, r) - y"{0, r) = iz?(0, r). (17) 

Nevertheless, when a is close to 1, namely a = 1 — Ce'' with r > and C > 0, the perturbative term 
in (|3]) has a non uniform bound for x € [0, 27r]. Indeed, its maximum (in absolute value), which is 
fie^ /{I — a)^ = /ie''~^''/C^, is reached at x = 37r/2. Therefore, it is natural to expect that the invariant 
manifolds of the perturbed system remain /xe''-close to the unperturbed separatrix only before they reach 
a neighborhood of the section x = 37r/2. For this reason, in this case we will measure the distance in this 
section, where one expects that the perturbed manifolds are closer. This section, by (J13p . corresponds to 
u = ln(l + v2)- Therefore, we define the d{T)^ the distance between the perturbed invariant manifolds 
at the section x = 37r/2, which is given by 

d(r) = y' fin (l + Vf] .t) - y" fin f 1 + V2) , r) = 2y/2D fin f 1 + \/2) , r) . (18) 



For a close to 1, one would expect that the distance d{T) is bigger than d{T) since the manifolds deviate 
from the homoclinic after they cross the section x ~ 3tt/2. Nevertheless, when a is not close to 1, both 
quantities are equivalent. 

2.1 The Melnikov function 

In this section we compute the Melnikov function associated to Hamiltonian ([S]) with /(r) = sinr and 
wc study its dependence on e and a. Let us recall that the Melnikov function (|15p is given by 

.,, , f^°° sinh(u + s) cosh(w + s) / s\ , , ^ 

7W(u,r;e,a) =4 / L___! L___! ^sm(T+-\ds. (19) 

•^-oo (cosh (u + s) — 2asinh(ii + s)) ^ 

As it is well known, the evaluation of this integral can be done using Residuum Theory. To this end, one 
has to look for the poles of the function 

/3(u) = Binh(^)eosh(u) ^ ^^0) 

(cosh^(?i) — 2asinh(ii)) 

closest to the real axis. For a = 0, /3(ii) has order 3 poles at p^ = ±i7r/2. Nevertheless, for a > 
these poles bifurcate into a combination of zeros and poles. Since the size of the Melnikov function P^ 
depends on the location of these poles, one has to study their dependence on a. We will see that the 
relative size between the parameter a and the period 27r£ of the perturbation will lead to a significantly 
different size of (fT9)) . 

Considering the denominator of /3(u), namely (cosh (u) — 2asinh(u))^, one can easily see that the 
poles of (PU)) are the solutions of 

sinhu = azbivl ^ ct^, 
This equation has four families of solutions given by p- + 27rfci, p_|- + 27rfci, which are solutions of 



sinhw = a + iyl — a^ . 

and their conjugate families p_ + ^-nki and p^ + l-nki for k G Z. The singularities p_, p+, p_ and p_|_ are 
the closest to the real axis (see Figure [2]) and they satisfy 

TT 

< Imp_ < — < Imp+ < TT, 

forae (0,1). 

In particular, if one considers a as a small parameter, they have the expansions 

P±(a) =i^±(-l + »)\/a + C'(a) as a ^- 0. (21) 

In the other limiting regime a — > 1, their expansions are 

p_(a) = ln(l + V2) + z(l - a)i/2 + 0(1 - «) 

^ , as a ^- 1. (22) 

p+{a) ^ - ln(l + V2) + TTi ~ i{l ~ af^ + 0(1 - a) 

Then, one can see that the function /3 for — tt < |Imu| < tt and — 1 < Reu < 1 behaves as 

n, ^ (u-i7r/2)(u + zV2) 

^^^^ (u^p^Y{u-p^Y{u--p_Y(u--p^Y- ^ ^ 

Using this fact, one can easily compute the size of the Melnikov function in ((T9)) using Residuum Theory. 
Proposition 2.1. There exists Eq > such that e G (0,eo) a^cf a E (0, 1), 



Imu = TT 



+ P4 



Imu = TT 



P++ 






+p- 



P++ 



Ini U = TT 



+«f 



+-» 



+ /7_ 



+p- 



+p 



±. j. 



Imu = TT 



Figure 2: Location of the singularities u ^ p^, p+,'p_,'p_. In the left picture we show them for a small, 
in such a way that they are close to the singularities of the separatrix of the pendulum u = zLiTr/2. In 
the right picture we show them for a close to 1. Then, p_ and p_ approach the real axis while p+ and 
p_l_ approach the lines Imu = tt and Imu = — tt respectively. 



• If a satisfies < a < Ce'^ where v > 1 and C > are constants independent of e, the Melnikov 
function (|19p satisfies the asymptotic formula 



47r ZL/ /a 2L\\ 

M{u,T;e,a) ^ —e 2s ( cos(r - u/e) + O f — , e 2s jj. 



(24) 



• If a = a+e^ + O (e"^) for some constant a^, > 0, the Melnikov function (J19p satisfies the asymptotic 
formula 

u/e) + ( -e^^ 



M {u,T;e,a) = —e 2e sin (r - 



(25) 



for certain constants A = A(a*) and 0* — 0*(e, a*). 



• If a satisfies Ce'^ < a < 1 where v € (0,2) and C > are constants independent of e, the Melnikov 
function ([T5| satisfies the asymptotic formula 



Ai {u,T] e,a) 
where Si{a) are given by 



(52(a) 



+ Si{a) 



Im p— 



£ sin (r - - u/e) + [e 2e 



Siia) 
62(0.) 



2tt 
(1 - a2)3/2 V 
27r sinh /9_ 
(1 — a'^) cosh/9_ 



I sinh/9_ —Ml 



2\l/2 



(26) 

(27) 
(28) 



with sinhp_ = a + iy/l — oP- , and = (/)(£, a). 
When a ^ e^, the first statement of Proposition 12.11 ensures that the Melnikov function is non- 



degenerate. In the case a ^ e^, one would have to analyze the behavior of the constant A(a*), which 
appears in formula (|25p . to check that it does not vanish. Even if this constant is computable in our 
example, we will not study it, since it is not the purpose of this paper to study this particular case. We 
just want to remark that for general a*, the first asymptotic order of the Melnikov function has changed 



from (P^ to (psj) . In fact, as can be seen in formula (j26p . a ^ e^ is the transition value in which the 
size of the Melnikov function changes drastically. Indeed, formula (j26p and the expansion of p_ given in 
(P^ show that the size of the Melnikov function becomes bigger when a approaches 1. This particular 
behavior will be analyzed in Section 12.1.11 

Next corollary, whose proof is straightforward, analyzes the behavior of formula (j26p for a in the 
range £^ ^ a < ao < 1 for any fixed ao G (0, 1). Let us observe that p- has the asymptotic expansion 
((2T|) when a -^ Q. Next corollary shows that the Melnikov function is exponentially small for this range 
of the parameter a. 

Corollary 2.2. Let us fix any ao S (Ojl)- Then, the Melnikov function in (|26p has the following 
asymptotic formulas for e^ ^ a < ao < 1 • 



// a satisfies a ~ Ce'^ for certain v G (0, 2) and C > 0, 

(^gl + i//2 



M{u,T-e,Ce-) = -^^-^e ^ (sin(T - - ^e) + O (e''/')) , 



where 

Therefore, in this case the Melnikov function is non-degenerate since (^2 7^ 0. 
• If a Cz (0,ao] is independent of e, 

M {u, r; e, a) = ^^^i^e — (sin (r - (/. - u/s) + O (e)) . 

e 

It can be checked that (^2(0) ^ /or any a G (0, 1), and therefore, in this case the Melnikov function 
is also non-degenerate. 

Remark 2.3. // one takes a ^ e^ and v ^ 2 in formula (|26p does not obtain formula (|25p . The reason 
is that formula (j26p only takes into account the residuum of p_ since the residuum of p+ is exponentially 
small with respect to the one of p^. Nevertheless, this is not the case when a ~ e^ (see (|2ip ). Then, 
the residuums of both singularities p_ and p_|_ make a contribution to the Melnikov function of the same 
exponentially small order 0(e~^e~'^'^^). 

Remark 2.4. In the first statement of Provosition [Kl[ one can see that, when the strip of analyticity 
^ is wide enough, namely taking a ^ e^ with v > 1, the Melnikov function at first order behaves as 
in the entire case a = lGel97a^ . In other words, the Melnikov function does not realize the finiteness 
of the strip of analyticity and the exponentially small coefficient is given by the imaginary part of the 
singularities of the separatrix. 

On the other hand, when the strip of analyticity ^ is independent of £ or not extremely big with respect 
to £, namely taking a ^ e" with v e [0,2), the exponentially small coefficient does not coincide with the 
imaginary part of the singularity of the unperturbed separatrix. Instead it is given by the imaginary part 
of this new singularity p_ , which appears when one evaluates the perturbation along the unperturbed 
separatrix. Note that even if one takes, for instance, a ~ e, the strip of analyticity is of order ln(l/£) 
and one has that 

M {u,T;e,e) ^ e~'^e 2e . 

That is, even for perturbations with a wide strip of analyticity with respect to e, it appears a correcting 
term in the exponential. 

The case a ^ e^ is the boundary between these two different behaviors. In this case, the exponential 
coefficient is given by the imaginary part of the singularity of the unperturbed separatrix but the constant 
in front of the exponential is not the one given in the first statement of Proposition \2.1\ but a different 
one, which is given in the second statement for the value a ^ a^e^ + 0{£^). 
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Proof of Proposition \2.1\ The classical way to compute the Melnikov function is to change the path of 
integration to Ims — ±7r and apply Residuum Theory, using that the function j3 has order two poles at 
p_, p+, p_ and p_|_ in the strip — tt < Ims < tt. Nevertheless, to prove the first statement, instead of 
computing the residuums of both p± and p_|_ and look for the cancellations, we just expand the Melnikov 
integral in power series of a, since it is uniformly convergent as a real integral. Thus, we obtain. 



Al(M,r;£,a)=4^(n + l)2"a" / 



sinh' [u + s) 
{u + s) 



cosh' 



2n+3 



.in (r + J) 



ds. 



Then, the term for ?i = gives the first asymptotic order and can be computed using again Residuum 
Theory. To bound the other terms goes as follows. We Fourier-expand the terms in the series of the 
Melnikov integral in r and we change the path of integration to Ims = 7r/2 — e or Ims = — (7r/2 — e) 
depending on the harmonic. This gives us the exponentially small term e~5?. To bound the other terms 
in the integral, we use that 



sinh"+'(u + s) 



cosh 



2n+l 



(w + s) 



<(f 



2n+l 



for certain constant K > Q independent of e. Moreover using that 

1 



cosh (w + s) 



decays exponentially as Re s —^ ±00 and that it has poles of order two u + s — ±i7r/2, one can see that 



1 



cosh'^(u + s± (i7r/2-£)) 



rds < —. 



Then, one has that 



+°° sinh"+Vu + s) / s, , 

^ ' sm ( r + - ) ds 



00 cosh2"+3(u + s) 



< 



2n+2 



2e . 



Therefore, the remainder can be easily bounded provided < a <C e^. 

When a ~ e^ we cannot ensure that the preceding series is convergent and therefore we use directly 
Residuum Theory in the strip < Ims < n or — tt < Ims < depending on the harmonic. For instance, 
for the positive harmonic, we take the strip < Ims < tt, which contains the singularities p_ and p+, 
that are at a distance of order 0{e) from i7r/2 (sec (pijl ). Then, it is enough to compute the residuum of 
/3(s)e'*/'^ at these singularities, which using (P^ satisfies 



Res(^(s)e''/",s = /9±) = 



A± _ 


TT 


-e 


TT 

~2e 


e^ 




\ e 





for certain computable constants A± g C independent of e. 

In the case a ^ e'^ with < z^ < 2, we have that 7r/2 — Imp„, Imp^ — 7r/2 - 



e"/"^ > £. Therefore, 



Im p+ TT Iinp- 



Using these facts, we can compute the Melnikov integral changing the path up to Im s = 7r/2 and we just 
need to consider the residuum at the singularity /)_ , which can be explicitly computed. D 

Remark 2.5. All the singularities of the function f3{u) in (|20l) have different imaginary part for any 
a G (0,1). This is one of the reasons of the choice of the perturbation in (jH]), since this is not always the 
case. Let us consider, for instance, the pendulum with a different perturbation as 



fie 



sin a; t 

■ sm - , 



(1 — acosx) 



(29) 
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which has Hamiltonian function 



H [ x,y,- \ — - — h cos a; — 1 + ^e^'— 7 -sin- . 

e J 2 a(l — Qfcosa;) e 



Then, the Melnikov function is given by 



.,, t\ ^ f^°° cosh(M + s) sinh(u + s) . ft + s 
7W u, - = 4 / ^ ^ '-^ sm ' 



s^ -l-oc ((1 -Q;)cosh^(u + s) + 2q;)' 

To compute its size, one has to study the singularities of 

~, , coshwsinhw 

I3{u) - -2- 

((1 — a) cosh u + 2a) 

The closest ones to the real axis with positive imaginary part are 



p± ^ t— ± arcsmni 



2 V 1-a 

and therefore both have the same imaginary part for any a G (0, 1). In these cases, the Melnikov function, 
which is given by the residuums of both singularities, has the same size for any a satisfying < a < ao < 1 
for any a^ independent of e, but is divergent for a~\. 

2.1.1 Narrow strip of analyticity: a drastic change on the size of the Melnikov function 

In Corollarv l2.2l we have seen that the Melnikov function is exponentially small with respect to e provided 
< a < ao < 1 for any fixed ao G (0, 1). We devote this section to study this function when a is close to 
1, which will lead to a Melnikov function which is either exponentially small with a different exponential 
dependence on e or even to a non-exponcntially small Melnikov function. 

We consider a = 1 — Ce'' with r > and C > 0. In this setting, one can see that the analyticity 
strip (O of the Hamiltonian system (JH) is very small, of order ©(e'"/^). As a consequence, the singularity 
p~{a) is very close to the real line. Indeed, from ([2^ . one has 

Imp_(a) = Ci/2e'^/2^0 (£'■). 

Corollary 2.6. If a = \ — Ce'' with r > and C > 0, the Melnikov function in (|26p has the following 
asymptotic formulas 

• // < r < 2 

M (u, r; £, 1 - Can = ^ ^^,^ e — (sin (r - - u/e) + O [e^/', e'-^^')) . 



^/2Cel+'' 
• Ifr = 2, 



-VC 



M {u, r; s, 1 - Ce') = -^^j^ (sin (r ~ <j> - u/e) + 0[e)) , 

• Ifr>2, 

Al(., r; s, 1 - CO = ^=^^^ (sin (r - ^ - ./e) + O (e'V-^) ) , 

where (j) = (j){e,a) is the constant given in Provosition \2. 11 

Therefore, in all these cases the Melnikov function is non-degenerate. Moreover, in the last two cases, 
it is non- exponentially small. 

Remark 2.7. To illustrate which size has the Melnikov function for the range of a considered in the first 
statement of Corollaru \2.b\ we can take a = 1 — £. Then, using the expansion (j22p . one can see that 

_ 1 
M{u,T;e,l- e) ^ e~'^e V^. 

Namely, the Melnikov function is still exponentially small but it has a different exponential dependence 
on e. 
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2.2 Validity of the Melnikov function 

Once we have computed the Mehiikov function in Proposition 12.11 provided < a < Q;o < 1 for a fixed 
ao, we can compute the prediction it gives for the distance between the manifolds at the section x = tt, 
which we call do{T). Recall that by ((T5)) it is the first order in /x of the function ^(t) given in (|17p . 
If a satisfies < a < Ce"^ where i/ > 2 and C > are constants independent of e, dQ{T) is given by 

do(T) = 27r/ie'' ^e 2e cost + O — ,e 2e 



If a = a*£^ + 0{e^) for any constant a* > 0, it is given by, 

doir) - ^^/^e^-^e-^ sin(T - 0*) + O Le'^'^e"^) . 

Finally, if a satisfies Ce'^ < a < ao < 1 where i^ e (0, 2), C > and ao E (0, 1) are constants independent 
of £, it is given by 

do{r) = ^^Y^iJe"'-^e ^ sin(r - 0) + O ( Aie"e — 

As we have already explained, direct application of Melnikov theory only ensures that the distance 
between the manifolds is given by 

d(T)=do(T)+0(A"')- 

Therefore, the Melnikov function is the first order of the splitting provided /i is exponentially small 
with respect to e. However, it is well known that often, even if Melnikov theory cannot be applied, 
the Melnikov function is the true first order of the splitting (see |HMS88[ IDS921 IGel941 IDS97[ IGel97al 
IBF04[[BF05llGelOO[IGOS10[IBFGSllO . Next theorem shows that, under certain conditions, the Melnikov 
function gives the true first order of the distance between the manifolds. We want to point out that the 
only available proof of the correct prediction of the Melnikov function for mcromorphic perturbations 
is jGel97a| . in which the author considers systems with analyticity strip wide enough with respect to e. 
All the other references deal with polynomial perturbations. Thus, this theorem is, as far as the authors 
know, the first one which shows the dependence of the size of the splitting on the width of the analyticity 
strip. 

Theorem 2.8. Let us consider any fiQ > and ao € (0,1). Then, there exists Eq > such that for 
|/i| < fiQ, e G (0,eo) 0,'iT'd a G (0,ao] such that e^^^{e + ^/a) is small enough, 

• If a satisfies < a < Ce'^ where v > 1 and C > are constants independent of e, the invariant 
manifolds split and their distance at the section x — "k is given by 

d{T) =2Ti^e'^-^e-Te (cost + O f-^^Me" 



• If a — a*e^ + 0(e'^) for some constant a* > and the constant A(a*) introduced in Provosition \2.1\ 
satisfies A(a, ) ^ 0, the invariant manifolds split and their distance at the section x ^ ir is given by 

d{T) ^ \^^^^^e-^-^e-fe (sin(r - 0*) + O (^e")) • 

• If a satisfies and Ce'^ < a < a^ where v G (0, 2] and C > are constants independent of e, the 
invariant manifolds split and their distance at the section x ~ t: is given by 

d{T) = M^M^e'f-ig-^ (sin(T - 0) + O (Ai£""iyS,e)) . 
The proof of this theorem is deferred to Section 2] 
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Remark 2.9. // one considers a ^ e^ with v > 2 the condition e''~^(e + -y/a) small enough is equivalent 
to requiring rj > 0. This condition is the same that has to be required if one considers the polynomial case 
a = (see \Tre9'7\ IBFGSll]/ ). In other words, in the case v > 2 the validity of the Melnikov function is 
the same as in the polynomial case. 

On the other hand, if one considers a ^ e'^ with v < 2, the condition e^^^{e + ^/a) small enough is 

r/ - 1 + ^ > 0. 

In particular, if one considers a as a parameter independent of e, one has to require 77 > 1. 

The limit cases in the previous settings, rj = and rj — 1 + v 12 = respectively, are usually called 
singular cases (see WFGSllf ). In Appendix\^ we make some remarks about how these cases could be 
studied. 

2.2.1 Narrow strip of analyticity: validity of the Melnikov function 

In Section 12.1.11 we have seen how the Melnikov function changes drastically its size when a = 1 — Ce'' 
with r > 0, even becoming non-exponentially small if r > 2. In this section we prove that it gives the 
correct first order for the splitting of separatrices. Note that, in the range of parameters for which the 
Melnikov function is not exponentially small, one can just apply classical perturbation techniques to 
prove that this function gives the first order of the splitting. 

First, we give the prediction of the distance given by the Melnikov function when a = 1 — Ce'' with 
r > and C > 0, which can be deduced from Corollarv l2.6l Recall that, as we have explained in Section 
[21 in this case we study the distance between the manifolds at the section a; = 37r/2. We have called ^(t) 
to this distance (see dH])) and we call do{T) to the Melnikov prediction of this distance. 

If r G (0, 2), calling u* = ln(l + V2),the Melnikov prediction is given by 

do{T) = ^£"— 1 — e — (sin (r - - u,/e) + O (^e''/^,s'-''/^j) , (30) 

whereas if r = 2 is given by 

doir) ^ ^£"-3_-^ (sin(r - - u*/e) + 0{e)) . (31) 

Finally, if r > 2 is given by 

doir) = ;ie"-3'-/2^ (sin(r -<l>-u*/e) + {e''/'-')) ■ (32) 

In these cases, direct application of Melnikov theory only ensures that the distance between the manifolds 
is given by 

d(T)-do(r)+o(MV(''-2-))^ 

since the perturbation has size 0(/i£''~^''). Next theorem widens the range of the validity of this prediction 
under certain hypotheses. 

Theorem 2.10. Let us consider any fxo > and let us assume rj > niax{r + 1, 3r/2}. Then, there exists 
Eo > such that for |/i| < ^q, e G (0, Eq) a'^'^ a = 1 — Ce*" with r > and C > 0, the invariant manifolds 
split and their distance on the section x = 37r/2 is given by. 



• If0<r<2, 

d(r) = ;ie"-'^-i^e"^ (sin(r - (f> - u*/e) + O (^^ie'^-''-\e^-'-^^,e''/^^^ 

• Ifr = 2, 

,2 

C3/2 



2Trp-VC 

^(^) = ^^^'"'^^37^ («in(T -cb- u*/e) + O (e, /ie"-^)) 
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• Ifr>2, 

d[r) = M^"-'''/' J^ (sin(r - - u*/e) + O (e^'^-\ iie^-^^/^)) , 

where u* = ln(l + \/2) and <f> ~ (j>{e, a) is the constant given in Corollary \2. 61 

3 The Melnikov function and its validity in the quasiperiodic 
case 

Analogously to the periodic case, we work in the extended phase space that now is {x, y, 61,62) G T x M x 
T^. Note that we do not work in the whole symplectic space, namely we omit the actions (/i,/2), since 
they do not have any dynamical interest. In this setting T = {(0, 0,6*1, 6*2); (^1,^2) S T^} is a normally 
hyperbolic torus and, for fJ. = 0, its stable and unstable invariant manifolds coincide along the homoclinic 
manifold 

W^iT) =W%T) ^ {{x,y,6,,62); Ho{x,y) ^ 0} 

= {{x,y,6i,62) = {xo{u),yo{u),6i,62y, (^,^1,^2) G R x T^} , 

where {xo{u), yo{u)) is the time parameterization of the homoclinic orbit, which is given in (|13p . 
Then, we look for the perturbed manifolds as 

J a- = xo{u) 
\y^y'^^%u,6,,62). 

Here (u, 61, 62) € (—00, [/) x T^, for certain [/ > 0, for the unstable manifold and (u, 6*1, 62) G {—U, +00) x 
T^ for the stable one. Again, the Lagrangian character of the manifolds implies that the functions y"''* 
are given by 

y^'^ (11,61,62) = -^duT"''{u,6i,62). 

Therefore, as in the periodic case, a natural way to measure the difference between the manifolds is to 
compute 

D{u,6i,62) ^ duT'iu,6i,62) ~ duT''{u,6i,62), 

whose first order in n is given by the Melnikov function 

M {u, 61,62) = / {Ho, Hi} i^xoiu + s), yo{u + s),6i + -, ^2 + 7- j ds. 

Namely, one has that 

Diu, 61,62) - fieUIiu, 61,62) + O {^l^e^'') . 

As in the periodic case, for a small or fixed independently of e, we measure the distance at the section 
X = TT, which is given by 

di6i, 62) ^y\0, 61,62) -y''{0,6i, 92)^^0(0,61,62), (33) 

whereas in the case a = 1 — Ce'' with r > and C > 0, we measure it at the section x = 37r/2, which is 
given by 

didi, 62) ^ y' (\n(l + V2) ,61,62) -y'' (in (1 + V2) ,61,62) = 2V2D (\n (l + V2) ,6^,62) . (34) 
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3.1 The Melnikov function 

If one applies the Poincare- Melnikov method to Hamiltonian P^ , obtains 

M{u,ei,92;e,a) = 4 f P{u + s)F (61 + -,9-2 + —) ds, (35) 

where /3 is the function defined in (^0)1 and F is the quasiperiodic perturbation ^ satisfying (fT(7|) and 

(ED- 

In the quasiperiodic case, it is a well known fact that the size of the Melnikov function is not given 
by its first harmonic (see |Sim941 IDGJS971 IDJSG99] ). Instead, the leading harmonic depends on e. For 
this reason, we need to compute carefully the size of all harmonics. We compute them using Residuum 
Theory and the properties of the function /3 given in Section[51 We follow the same approach of |DGJS97] . 
For this reason, let us first introduce certain functions and constants defined in that paper. 

We define the 2 In 7-periodic function c{6) defined by 

c{S) = Co cosh ( ^-^ ) for 6e[So~ \nj, 6a + In 7], (36) 

where 



C, = 2,m±^,So^lne* ^* - ±±J:^ 



7 + 7 ^ ' ' 7^(?'i7 + ?'2) 

and continued by 2 In 7-periodicity onto the whole real axis. Note that the constant Cq is slightly different 
from the one considered in [DGJS97| since in that paper, it includes a 7r/2 coefficient coming from the 
imaginary part of the singularity of the unperturbed separatarix. Since in this paper the singularity 
changes with respect to a, we have defined a new constant which is independent of it. Then, the formulas 
of the splitting, given in Proposition 13.11 and Theorems 13.51 and 13.71 will contain the dependence on the 
imaginary part of the singularity explicitly. Following the lines in JDGJS97] , one can see that the function 
c{6) oscillates and has lower and upper bounds independent of e. 

One can use the function c{S) to give the size of the Melnikov function as done in |DGJS97] . In the 
next proposition we see how this size changes depending on the relation between e and a. In all the 
results in the quasiperiodic case we include the case a = 0, since the proof we present is also valid in this 
case and slightly improves the results in the literature JDGJS971 ISauOl] . 

Proposition 3.1. Let us assume (|10p and (|lip and fix uq > 0. Then, there exists Eq > such that for 
e £ (0, Eq) and a e [0, 1), 

// a satisfies < a < Cs'^ where v > I and C > are constants independent of e, the Melnikov 
function (|35p satisfies that 

Ci ~c{ln{2eM).f§ M / C'2 -c(ln(2e/^))./^ 

— e V ^^ < sup \M{u,9i,62]e,a)\ < — e V ^^ (37) 

£ (u,ei,e2)6(-uo,"o)xT2 £ 

for certain constants uq,Ci,C2 > 0. 



• 



• 



In the intermediate case a = a*e + 0{e ), one can only give upper hounds of type 

C2 -c(lii(2e/7r))i/i|^ 

sup |A^ (m, 6'i,6'2;£, Qf)| < — e V ^^ 



(«,6»i,e2)6(-«o,i»o)xT2 

for certain constant uq,C2 > 0. 



£ 
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• If a satisfies Ce'^ < a < 1 — Ce'' where i' G (0, 1], r e [0,2) and C > are constants independent 
of e, the Melnikov function satisfies 



Cr 



-c(ln(e/Im p_ )) 



ea(l - a)5/4 



Im p- 



< sup |A^ (u,6'i,6'2;£,a)| 



< 



Co 



-c(ln(e/Imp_)) 



^(1 - a)5/4 



Iin p - 



(38) 



for certain constants Ci,C2 > 0. 



Proof. As in the periodic case, to prove the first statement, instead of computing the residuums of both 
p± and look for cancellations, we just expand the Melnikov integral, which is uniformly convergent as a 
real integral, in power series of a as 



where 



Mniu,ei,92;e)^4in + 1)2 



M {u, 61,62] e, a) = ^ a^Mn iu, 61,62; e) , 

T1=0 



(39) 



cosh' 



2?i+3 



(u + s) 



FI61 



-,02 + — Us. 

e e 



The function AIq was computed in (DGJS97] . and in that paper it was shown that it can be bounded 
as formula ([57]) . The rest of the functions A^„ can be bounded as follows. We Fourier-expand A^„ in 
(61,62) and we change the path of integration to Ims = tt/2 — ^/e or Ims ~ —(jij^ — ^/e) depending on 
the sign oi k ■ uj = k\ -\- 7/02. Then, we can bound each harmonic as 



M^^\u;e) 



< 



n+l 



F^h] 



\k-Uj\ f TT 



-v^ 



Therefore, proceeding as in |DGJS97] . one can see that A^„ {u, 61, 62; s) can be bounded as 

\Mniu,6i,62;e)\< 



A'\"+^ -c(ln(2eA))^/^ 



where the function c is given in (pS)) . 

Therefore, if e ^ a, the series (j39|) is decreasing, and therefore, the leading term is given by the first 
order. 

For the second and third statements, we just apply the Residuum Theory to the Fourier harmonics 
of the Melnikov function ((35|) . In can be easily seen that they have size 



7WW(m) 



j^W 



k ■ LJ 



52{a) + (Si (a) 



|/i;-Lj|Im p- 



(40) 



where 5i{a) and 52{a) are the functions defined in ([27]) and (pS)) respectively. For the the range of a we 
are considering, the leading term is 82, which satisfies 

1 
(52(a) 



y/a{l - a)' 



Finally one has to proceed as in |DGJS97] . Even if the method in that paper docs not apply directly 
since now p_ depends on e, it is enough to consider as a new parameter q = e/Imp-, which is still small 
since we are assuming Ce'^ < a <\~ Ce^ where i' G (0, 1], r S [0, 2) and C > and therefore Imp_ ^ e. 
With this new parameter q, it is straightforward to obtain the size of the Melnikov function with the 
techniques in JDGJS97) . D 
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Remark 3.2. Note that, as in the periodic case, the size of the Melnikov function depends strongly on 
a. When the strip of analyticity (|6]) is wide enough, that is a ^ e'^ with v > 1, the Melnikov function 
behaves as in the entire case a = jDGJSQT]/ . In other words, the Melnikov function does not notice the 
finiteness of the strip of analyticity and the exponentially small coefficient is given by the imaginary part 
of the singularities of the separatrix. Nevertheless, note that the condition is different. In the periodic 
case was needed a <g; e^ instead of a-^ e. 

When the strip of analyticity is independent of e or not extremely big with respect to e, namely taking 
a ^ e'^ with v S [0, 1), the exponentially small coefficient multiplying the periodic function c{S) does not 
coincide with the imaginary part of the singularity of the unperturbed separatrix. Instead, it is given by 
the imaginary part of this new singularity p_ , which appears when one evaluates the perturbation along 
the unperturbed separatrix. Note that even if one takes, for instance, a = ^/e, which gives a strip of 
analyticity of order 2 l^ ~ ■ one has that, using (j2ip , 



M {u,9i,92;e,y/e) ^ e 



^3 -c(ln(e/Imp_))./^(l-4^) 



That is, even for perturbations with a wide strip of analyticity with respect to £, there appears a correcting 
term in the exponential. 

Notice that the case a ^ e is the boundary between these two different behaviors. 

3.1.1 Narrow strip of analyticity: a drastic change on the size of the Melnikov function 

In this section, wc study how the Melnikov increases when the analyticity strip shrinks, namely when 
a ~ 1 — Ce'' with r > and C > 0. Next corollary gives upper and lower bounds for it in this case. 

Corollary 3.3. Let us assume (|10p and (jlip and fix uq > 0. Then, if one takes a = 1 — Ce^ with C > 
and r > 0, the Melnikov function in (j38p has the following upper and lower bounds 

• If0<r<2, 

Cl — c(ln(e/Imp_))y — \txr n n M 

r- i/2+5W4 ^ ^^V . < sup \Miu,0i,e2;e,a)\ 



. -;=j^e-*<-'-»->'v'^ («) 



• Ifr>2, 

'^^ < sup \Miu,ei,e2;e,l-Ce'')\< ^^ 

(u.ei.e2)e{-uo,uo)xT 



£3rV2 -^ ,...^.^,:rK.^..^.^r,2'' ' '^'^' ' ^'^ ^3r/2 



for certain constants Ci, C2 > 0. 

The first statement of this corollary is just a rewriting of the third statement of Proposition 13. II The 
second one is a direct consequence of formula (^0]) , if one takes into account the definition of Si and S2 
in ([271) and ^^ and the asymptotics for p_ in ((22|) . 



Remark 3.4. When the strip of analyticity shrinks, which ocurrs when a approaches 1, the Melnikov 
function increases exponentially. For instance, if one takes a = 1 ~ e, using (I22p . 

.. c(ln(E/Imp_)) 

M{u,6i,62]e,l-e)r^ e—^e e^'^ 

The exponentially small Melnikov function keeps increasing until the limiting case a ^ 1 — £^ when it 
becomes non- exponentially small. 
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3.2 Validity of the Melnikov function 

Once we have computed the Melnikov function in Proposition 13. 1[ wc compute the prediction it gives for 
the distance between the manifolds at the section a; = tt for < a < ao < 1 for any fixed ao 

If a satisfies < a < Ce" where v > \ and C > are constants independent of e, we can ensure that 
do satisfies 



Ci\ix\e'^ e V 2e < max |do(^i, ^2)! < C2|/i|e^ e V ^^ (42) 

If a = a^^e + C(e^) for any constant a^ > 0, do satisfies 

max Mo(0i,02)| < C2|//|£''-ie"'^'"^''^"^^V ^^ . (43) 

(ei,e2)eT2 

Finally, if Ce"^ < a < ao < 1 where i/ G (0, 2) and C > arc constants independent of e, do satisfies 



]= e V e < max 00(^1,^2) 

V" (ei,e2)GT2' 



^ ^^^^^_,(ln(e/Imp-))V^^ (44) 



In all three formulas Ci, C2 > are constants independent of a and e. 

As we have already explained, direct application of Melnikov theory only ensures that the distance 
between the invariant manifolds on the section a; = tt is given by 

d{9ue2) ^ do{0i,e2) + O (fi^e^^) . 

Therefore, the Melnikov function psp . is, in principle, the first order provided fi is exponentially small 
with respect to s. Next theorem shows for which range of parameters a, e and ry, the Melnikov function 
gives the true first order. In the quasipcriodic case, it is not known whether this range is the optimal one 
for which the Melnikov function predicts the splitting correctly. 

Theorem 3.5. Let us consider any ^q > and ao G (0,1). Then, there exists £0 > such that for 
\fi\ < fiQ, e G (0,eo) o.'iT'd a G [0,ao] such that e^^^^{y/e + y/a) is small enough, 

• If a satisfies < a < Cs'^ where v > 1 and C > are constants independent of e, the distance 
between the invariant manifolds on the section x = tt is given by 



d{e,,02)^do{e,,e2)+O L|2^2,-2^-(ln(e/Imp_))/ 



Iin p_ 



// a satisfies Ce'^ < o; < Q^o < 1 where v G [0, 1) and C > are constants independent of e, the 
distance between the invariant manifolds on the section x ^ tt is given by 



d{e,,e,)^do{0,,e,)+0 \ |^|2,2,-2^-(ln(e/Imp_))V^^ 



The proof of this theorem is deferred to Section [S] 

Remark 3.6. Note that even if this theorem holds true provided £'^^^{\/e + ^/a) is small enough, the 
Melnikov prediction is the true first order only if one imposes a more restrictive condition. 
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Comparing the size of the Melnikov prediction (j42p and the size of the remainder, in the case a <^ e, 
one has to impose the condition 77 > 1. Let us observe that when a = 0, this theorem, which gives that 
Melnikov provides the true first order of the distance if rj > 1, is an improvement with respect to the 
sharpest previous result jDGJSQT]/ , which needed the condition 77 > 3 (see also JSauOlf ). 

In the case a 3> £ one has to impose £^"^^[01 <C 1. 

Finally, note that the intermediate case a = a*£ + O(e^) is included in the first statement of the 
theorem. Nevertheless, if one would check that do is nondegenerate, Melinkov theory would give also in 
this case the correct prediction. 

3.2.1 Narrow strip of analyticity: validity of the Melnikov function 

In Section [3TO] we have studied the size of the Melnikov function when a = 1 — Ce^ with C,r > 0. We 
devote this section to show for which range of parameters the Melnikov function gives the correct first 
order of the splitting. 

First, wc give the prediction of the distance given by the Melnikov function if a = 1 — Ce'' with r > 
at the section x ~ 37r/2. If < r < 2, 



-c(ln(e/Imp_)))'^^^^^^ 



(ei,e2)eT2 



whereas if r > 2 



< C2|Mk''-^-'fe-=^'"'^'/'"''-»V^ (45) 



Ci|/i|£''""'^'< max doi0i,02) < CslMk""'"^ ■ (46) 

(ei.e2)eT2 

Next theorem proves the validity of this prediction under certain hypotheses. 

Theorem 3.7. Let us consider any /io > and a = 1 — Ce'' with r > 0, C > 0, and let us assume 
rj > max{r + l,3r/2}. Then, there exists eo > such that for \^\ < fj.Q and e G (0,eo), the invariant 
manifolds split and the distance between them on the section x = 37r/2 is given by, 



• // < r < 2, 



d(0i,02) = do(^i,e2) + ( |/.p£2"'-2'-2e-^'^'"^^/''"''-»V 



Im p- 



• Ifr>2, 

d{9i,e2) = do{ei,e2) + o {fie^^-'-') . 
4 The periodic case: Proof of Theorem 12.81 and 12.101 

In Section im we prove Theorem l2.8l that is when the parameter a is bounded away from 1. Nevertheless, 
we prove the result for a wider range of the parameter a. Therefore, our proof deals also with the first 
statement of Theorem 12.101 For this reason, during Section 14.11 we assume the condition 

< a < 1 - e^. (47) 

Note, that this range of parameters correspond to an exponentially small Melnikov function. The second 
statement of Theorem 12.101 is proved in Section 14.21 
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4. 1 Proof of Theorem [278] 

To study the splitting of scparatriccs wc follow the approach proposed in [LMS03[ ISauOll IGUSIO] , which 
was inspired by Poincare. Namely we use the fact that the invariant manifolds are Lagrangian graphs. 
This allows us to look for parameterizations of the invariant manifolds as graphs of the gradient of certain 
generating functions. 

First, let us point out that we will consider r as a complex variable to take advantage of the analyticity 
of the Hamiltonian with respect to this variable. To this end, we consider a fixed a > and we take 

reT^ = {reC/Z: |ImT| <<j}. 

Note that since the Hamiltonian is entire in r, we can take any a. From now on, all the constants 
appearing in this section will depend on a. As we will see during the proof, in the periodic case, the 
analyticity strip with respect to r does not play any role in the size of the splitting. 
Let us consider consider the symplectic change of variables 

X — xq{u) — 4arctan(e") 

w coshu (48) 

y = — r^ = — ^ — w, 

which was introduced in IBal06j (see also |LMS03[ ISauOl] ) and reparameterize time as r = t/e. With 
these new variables, the Hamiltonian function ([5]) reads 

sH{u,w,t) =eH ixoiu), — r^,T] (49) 

and the unperturbed separatrix can be parameterized as a graph as 

w = yQ{u) = — -3—. 
cosli u 

Then, one can look for the perturbed invariant manifolds as graphs of the gradient of generating functions 
T"'*(u,t). Namely, we look for functions T"'''(u,r) such that the invariant manifolds are given by 
w = duT^'^iuTT). Moreover, these functions are solutions of the Hamilton- Jacobi equation 



This equation reads 



with 



e-^a^T + H (w, a„r, r) = 0. (50) 



e-^drT + £^!LJ^ (a„T)' ^ + iis'^^iu) sinr = (51) 

8 cosh u 



*(u)=V(xo(«)), (52) 

where i\} is the function considered in ([S]) and xq(u) is the first component of the time-parameterization 



of the separatrix (see 

Moreover, we impose the asymptotic conditions 

lim y^^{u) ■ duT^iu,T) = (for the unstable manifold) (53) 

Rg u— f — 00 

lim y'^^iu) ■ duT^{u,T) = (for the stable manifold). (54) 

Re li— f +00 

One can easily see that for /^ = 0, the solution of equation ([FT]) is just 

e" 
To(u)=4^— , (55) 

coshw 
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which is the generating function of the unperturbed separatrix and has singularities at u = ±i7r/2. 
Nevertheless, as we have explained in Section [21 the term ^(u) has singularities ai u = p±,'pj.. In 
particular, u = p-,'p_ are closer to the real axis than u = ±i7r/2, u = p+ and u — p_|_. Therefore, 
to study the exponentially small splitting of separatrices for this kind of systems, one has to look for 
parameterizations of the invariant manifolds, namely solutions of equation (j51|) . in complex domains up 
to a distance of order 0{e) of the singularities u = p-^'p_. To this end, we define the domains 



13" = {u e C; jlmMJ < — tan/3i (Reu — Rcp_) + Im/5_ — ke} 
D^. = {u E C; |Imu| < tan/3i (Reu — Rcp_) + Imp_ — ke} , 

where /3i > is an angle independent of e and k > is such that Im/9_ — ne > (see Figure [3]). 



(56) 





Figure 3: The domains D" and Df^ defined in (j56p . 

Note that this is a significant difference in comparison with the papers dealing with the exponentially 
small splitting of separatrices problem for algebraic or trigonometric polynomial perturbations (see |DS971 
IGel97bl IBF041 IBFGSll] ). In these cases, one needs to look for the parameterizations up to a distance 
of order 0{e) of the singularities of the unperturbed separatrix. Nevertheless, in this paper, since the 
perturbation changes the location of the singularities of the perturbed invariant manifolds so drastically, 
one has to study the invariant manifolds close to these new singularities. 

Let us recall that when a is close to 1, namely, when the analyticity strip of the Hamiltonian is very 
narrow, the singularity p_ is given by 

p_ = ln(l + \/2) + i{l - a)i/2 + 0(1 _ a). 

Therefore, Imp_ <^ e ii 1 — a -^ e^. In this case, as the Melnikov function is not exponentially small 
(see Corollarv l2.6p . we can use a classical perturbative approach to prove its validity. We leave this easier 
yet surprising case to Section 14.21 Thus, from now on, we assume the condition (j47p and we proceed to 
prove Theorem 12.81 This condition ensures that the intersection D]^ n D'^k n M contains a fundamental 



domain, since it is of size 0{\/l — a). Note that if a is bounded away from 1, one can choose the angle 
/3i so that the domain D^ n D'^k D M contains the point u = (which corresponds to the section x = tt). 
If a = 1 - Ce'' with r e (0, 2), the domain D]^ n D^k n M contains the point u = ln(l + V2) (which 
corresponds to the section x = 3tt/2). One could change slightly the domain so that it would contain 
u = also in this latter case. Nevertheless, in this case one would have to take j3i ~ e''/^. This would 
lead to worse estimates which would require a stronger condition on 77 for the Melnikov function to be 
the true first order of the splitting. 

The next theorem gives the existence of the invariant manifolds in the domains D* with * = u, s 
defined in (|55)l . We state the results for the unstable invariant manifold. The stable one has analogous 
properties. 
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Theorem 4.1. Let us fix ki > 0. Then, there exists eg > such that for e e (0,eo), a £ (0,1), 
fj, e B{iiq), if ^i-z~^e''~^ is small enough and (|47p is satisfied, the Hamilton- Jacobi equation (j50p has a 
unique (modulo an additive constant) real-analytic solution in D" xTo- satisfying the asymptotic condition 

dSSD. 

Moreover, there exists a real constant bi > independent of e, /i and a, such that for {u,t) e D"^ xTg-, 

fell/lie"-! 



\duT"^iu,T)-duToiu)\< 



(£ + y^)(l-a) 



\dlT-{u,r)^dlUu)\< ^^l'^"''"' 



(e + Va)(l-a)' 
Furthermore, if we define the half Melnikov function 

,.v/ ^ . [^ sinh(u + s) cosh(u + s) . / s\ , ,^„, 

M {u,t) ^-A ^ ^ ^smlT+ -jds, (57) 

J-oo (cosh (m + s) — 2asinh(u + s)) ^ 

the generating function T" satisfies that, for {u,t) G Z?"^ x T^-, 

2r,-2 

|9„T"(u,T)-a„ro(")-M^''-M"(",T)| <6i|Aip- r^. (58) 

(1 — ay 

The proof of this theorem is deferred to Section 14.1.11 The parameterization of the stable manifold 
has analogous properties. In particular, we can define 

,_, , ^ f^°° sinh(it + s) cosh(w + s) . / s\ , ,^„, 

M''iu,T) =A ^ i ^sm[T+ -\ds, (59) 

"'o (cosh^(M + s) - 2asinh(u + s)) ^ ^^ 

and then, for (w, r) e -D^^ x T^-, 

2ri-2 

\duT'iu,T) - duToiu) - fie'^M-^u^T)] < \^i\^- -. (60) 

(l-a)2 

Note that the Melnikov function defined in ([T^ is simply 

M{u,T;e,a) = A^*(u,r) - A<"(u,r). 

From now on, we omit the dependence on e and a of the Melnikov function M, which we denote by 
M{u,t). 

Once we know the existence of parameterizations of the invariant manifolds, the next step is to study 
their difference. To this end, we define 

A(u,r) =r''(u,T)-r"(u,T). (61) 

This function is defined in R^^ x To-, where i?^, for any k > 0, is the romboidal domain 

R. = D:n Dl (62) 

Subtracting equation (|?T|) for both T'^ and T", one can easily see that A G Ker£j for 

C, = e-^dr + (^2^ {duT'{u,T) + duT^{u,T))\ a„. (63) 

Since Theorem 14. II ensures that the perturbed invariant manifolds are well approximated by the unper- 
turbed separatrix in the domains D"^ and Z)^^ , we know that the operator C^ is close to the constant 
coefficients operator 

C,^e-^dr + du. (64) 

As it is well known, any function which is defined for (u,r) G {u G C : Reu = a,Imu G [— ro,ro]} x T„ 
and belongs to the kernel of £e. it is defined in the strip {|ImM| < rg} x T^ and has exponentially small 
bounds for real values of the variables. This fact is summarized in the next lemma, whose proof follows 
the same lines as the one of the slightly different Lemma 3.10 of |GOS10| . 
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Lemma 4.2. Let us consider a function Q{u^ r) analytic in {u, r) e {u e C : Re u ~ a, Im u G (— ''Oj ''o)} ^ 
To- which is solution of C^C = 0. Then, C, can he extended analytically to {|Iinw| < tq} x To- and its mean 
value 



1 c27r 



(c> = ^y^ au,r)dT 

does not depend on u. Moreover, for r g (0, rg) and a' € (0, cr), we define 

Mr— max _ \duC{u,T)\ . (65) 

(u,r)G [~^r,^r] xT^/ 

Then, provided e is small enough, for (u, t) G M x T^-' t/ie following inequality holds 

\duCiu,T)\ <4M,e-f. 

To apply this lemma to study the difference between the invariant manifolds, we follow |Sau01| (see 
also |GOS10| ). Namely, we look for a change of variables which conjugates £e in (|55|) with C^ in (|M)) . 

Theorem 4.3. Lei us consider the constant ki defined in Theorem\4-l\ and let us fix any k^ > K2 > ni- 



-1 e+y/g • 



1-Q 



JS 



Then, there exists Eq > smc/i i/iai /or e e (0,£o) a'^'^ Q^ G (0, 1) satisfying (|47p anrf f/iaf e'' 
small enough, there exists a real- analytic function C defined in i?^^ x T^- such that the change 

{u,t)^{v+C{v,t),t) (66) 

conjugates the operators C^ and C^ defined in (|63p and (j64[) respectively. Moreover, for {v, r) 6 i^^a x To-, 
V -\-C{v,t) € i?Ki OLnd there exists a constant 62 > such that 

\ — a 

1 — a 

Furthermore, (u,t) ~ {v+C{v,t),t) is invertible and its inverse is of the form {v,t) = {u + V{u,t),t) 
where V is a function defined for {u,t) G Rk^ x T^, which satisfies 

1 — a 
and that u + V(u, r) G i?^^ for [u, r) G Rk^ x To-. 

The proof of this theorem is deferred to Section 14.1.21 

Once we have obtained the change of variables C, we are able to prove the validity of the Melnikov 
prediction. We prove it by bounding 9„A(M,r) — Ai{u,T), where A and A^ arc the functions defined in 
([CT|) and ([T^ respectively. 

As a first step, we consider 

d,{^{v+C{v,T),T))~M{v,T) 

where C is the function obtained in Theorem 14.31 

Theorem 4.4. There exists Eq > and 63 > such that for any e G (0,eo) and a G (0, 1) satisfying 
(j47p and that e^'^^ ^ i-a *'^ small enough, the following hound is satisfied 



\d,{^{v+C{v,T),T))-M{v,T)\<h^\tJL 



2 



(l-a)2 

for V G i?K2 n M and t G T, 
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Proof. First, we define the Melnikov potential L (see |DG00j ). namely a function such that dvL{v,T) 
Ai{v, r), and we study A{v +C{v, r), t) — L{v, r). 

As A e Ker££, where C^ is the operator in (j63p . by Theorem 14.31 the function $(w,t) = A{v 
C{v, r), t)~L{v, t) S Ker/Cg, where £e is the operator defined in (j64[) . Therefore we can apply Lemma[ 
To this end, we have to bound 9„$(u, r) in the domain R^^^ xlTcr- We split $ as $(u, t) = ^i{v, r)+$2(^', t) 
where 

$i(u, r) = A (w + C(u, r), r) - A(u, r) 
$2(u, r) =: A(u, r) - L(w, r). 

where A is the functions defined in (pT|) and L is the Melnikov potential. 

To bound 9t,$i, one has to take into account that A = (T^ — Tq) — (T" — Tg) and therefore, it is 
enough to consider the bounds obtained in Theorems 14.11 and 14.31 to obtain 

\d.<i>iiv,T)\<K\fi\e^''-'- ^ 



(1 - a)2 

for (w,t) € Rk,2 X TTcr. For the second term, it is enough to use bounds (|58p and (IBU]) to obtain 

1 



|a.$2(^,r)|<iC|M|e 



2r)-2_ 



(1 - a)2 



for {v,t) G i?„, xT^. 

Therefore, we have that 



\d,Hv,T)\<K\^i\e'''-^- ^ 



(l-a)2 
for (w, r) e i?K2 X Tg. and then, it is enough to apply Lemma [4.21 to finish the proof of Theorem 14.41 D 



From this result and the exponential smallness of M , given in Proposition 12. 1[ and considering the 
inverse change (v,t) = (u + V(u,t),t) obtained in Theorem 14.31 it is straightforward to obtain an 
exponentially small bound for duA.{u,T) — A4{u,t). It is stated in the next corollary, whose proof is 
straightforward . 

Corollary 4.5. There exists eo > and 64 > such that for any e G (0,eo) cii^d, a G (0,1) satisfying 
(|47p and that g'?-! ^^^^ is small enough, the following bound is satisfied 



£277—2 Imp_ 

\duAiu,T) - Miu,T)\ <bi\ti\''- -e — 

(1 — a)'' 



for u G i?K3 n M and t G T. 



From this corollary and using that, by Proposition 12.11 and Corollarv 12.21 the Melnikov function has 
non-dcgencratc zeros, one can see that if a is bounded away from 1 and satisfies that £''~^(e + y/a) 
is small enough, the manifolds intersect transversally and their distance satisfy the desired asymptotic 
formula. This finishes the proof of Theorem H^ (see jBFGSllj V 

4.1.1 The invariant manifolds: Proof of Theorem 14.11 

Since the proof for both invariant manifolds is analogous, we only deal with the unstable case. We look 
for a solution of equation ([5T|) satisfying the asymptotic condition (|53p . We look for it as a perturbation 
of the unperturbed separatrix Tq in (|55p and therefore we work with 

Q{u,t)=T{u,t)^To{u). (67) 

Replacing T in equation (|5ip . it is straightforward to see that the equation for Q reads 

C,Q^J^iduQ,u,T) (68) 
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where £e is the operator defined in (|M| and 



cosh^ 



J"(/i,u,t) = /i^-^£''*(u)sinr. (69) 

8 

where ^(w) is the function defined in ([5^. 

We devote the rest of the section to obtain a solution of equation (|55)) which is defined in _D" x T^. 
and satisfies the asymptotic condition ((55)) . 

We start defining a norm for functions defined in the domain _D" with k > 0. Since we want to capture 
their behavior both as Reu — >■ — oo, namely its exponential decay, and for u close to m = ia, we consider 
weighted norms with different weights. For this reason we consider fj > to divide D" by the vertical 
line Rew = —U. Then, given k > and an analytic function h : D]^ — > C, we consider 



||/i|| = sup |e-^"/i(u)| 

ueDi^n{Rcu<-u} 

+ sup \{u- p^)^(u- p+){u--p_f{u--p^)h{u)\. 

ueD^n{Rc u>-u} 

For analytic functions h : D" x T^ — > C, we consider the corresponding Fourier norm 



(70) 



\hL^J2h^'^ 



elfck. 



fcGZ 

and the following Banach space 

£^^ = {h : Dl X T^ ^ C; real-analytic, ||/i||„ < oo} . (71) 

To obtain the solutions of equation (|68)) . we need to solve an equation of the form C^h = g, where Cg 
is the differential operator defined in (|64p . Note that C^ is invertible in £k,o-- It turns out that its inverse 
is Q^ defined by 

.0 

g^{h){u,T)^ h{u + s,T + s-\s)ds. (72) 



Lemma 4.6. The operator Qg in (j72p satisfies the following properties. 

1. Q^ is linear from f^^o- to itself commutes with du and satisfies C^ o Ge ^ Id- 

2. If h e £K.,a, then 



Furthermore, if (h) ~ 0, then 
3. If h e £n,t7, then duQe{h) e ^k.o- and 



\mh)t < K\\h\\„. 

\geih)t < Ke\\h\\,. 



WduGeWt < K\\h\\,. 
Proof. It follows the same lines as the proof of Lemma 5.5 in [GOSIO] . D 



Once we have obtained an inverse of the operator ££ defined in (|64p we can obtain solutions of equation 
(j68p using a fixed point argument. Then, Theorem 14. II is a straightforward consequence of the following 
proposition. 

Proposition 4.7. Let us fix ki > 0. There exists Eq > such that for any e £ (0,eo) and a e (0, 1) 
satisfying (|47l) and that ^i_^ £^^^ is small enough, there exists a function Q defined in D^^ x T^ such 
that duQ G Stii.a is a fixed point of the operator 

T'i.h) ^ duGeJ'ih), (73) 

26 



where Q^ and J- are the operators defined in (|72p and (|69p respectively. Furthermore, there exists a 
constant bi > such that, 

\\duQL<bMe''^\ 

Moreover, if we consider the half Melnikov function defined in (|57p , we have that 



2 2r,f +\/" 



|a„g-/ie"A^"t<ii>|V''- 



(74) 



Proof. Let us consider kq < ki- It is straightforward to see that T is weh defined from Skq-u to itself. 
We arc going to prove that there exists a constant 6i > such that J^ sends B{bi\fj.\e^^^) C £^0.0- to 
itself and is contractive there. 

Let us first consider T (0). From the definition of T in ([75)1 . the definition of T in (|M)l and using 
Lemma [1751 we have that 



J-"(0)(w,t) = 9„ge-^(0)(u,r) = ^fie^g, {^'{u)smT) . 



(75) 



To bound it, first let us point out that \E''(u) = I3{u) where /?(u) is the function defined in (PO]) . We 
bound each term for Imu > 0, the other case is analogous. Using (P^ and taking into account that 



u — i7r/2 



u- P+ 



< 1, 



(76) 



one can easily see that ||^'(it)|| < K. Then, taking into account that sinr has zero average and applying 
Lemma [4.61 there exists a constant 61 > such that 



-^ (0) 



< ylA^k^^^ 



To bound the Lipschitz constant, let us consider /ii,/i2 G B{bi\^\e^^^) G Skq^ct- To bound ||J" (/i2) — 
•7^ {hi)\\cr, we need first the following bounds for u G D^o^ 



u — z7r/2 



u — ip- 



and 



<1 + 



1 



iTT/2 — p- 



u — ip- 



{u — ip^){u — ip^) 



< 



<1 + K 



K 

e(l — a) ' 



(77) 
(78) 



which are a direct consequence of (|21[) and (22). 
Then, it is easy to see that 



r{h2)-r{h{, 



< n^ ^ f 1 + — ) \\h2 + hiU\h2 - h,\ 
a e(l — a) \ £ 



I — a 



hlWa. 



Then, using that by hypothesis ^^_^ e'' ^ is small enough. 



LipJ" = /YHi±^e"-i<l/2 
1 — a 

and therefore J^ is contractive from the ball B{bi\p\£^'^^) C ^KQ.ff into itself, and it has a unique fixed 
point h* . Moreover, since it has exponential decay as Reu ^ — cx), we can take 

Q{u,t) = / h*{v,T)dv. 
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To obtain the bound for d^Q, it is enough to apply Cauchy estimates to the nested domains Z?"^ C -D"^ 
(see, for instance, {GOSlOj V and rename &i if necessary. 

Finally, to prove ([71)) . it is enough to point out that F (0) = fie'^Ai^ and consider the obtained bound 
for the Lipschitz constant. D 

4.1.2 Straightening the operator ^gt proof of Theorem 14.31 

We devote this section to prove Theorem 14.31 It is a well known fact, see for instance Lemma 6.3 
of |GQS10| . that looking for a change (pS)) which conjugates £e and C^ defined in ([55]) and (|M)) is 
equivalent to looking for a function C solution of the equation 

cosh u 



C,C{v,t) = ^^(a„T^(u,r) + 9„T"(u,T)) 



- 1. 

u—v-{-C{v,t) 

Taking into account the definition of Tq and Q in (|55p and (j67p . this equation can be written as 

C,C = J{C) (79) 

where 

cosh u 



J{h){v,T) = -^— {d^Q-{u,T) + duQ'^iu.T)) 



(80) 

ti— i;+/i(f ,t) 

To look for a solution of this equation, we start by defining some norms and Banach spaces. Given n G N 
and a function h : R^ — > C, we define 

\\h\U= sup \{v-p^T{v-Ti_rh{u)\. (81) 

Moreover for analytic functions h : i?^ x Tg- — > C, we define the corresponding Fourier norm 



l/^lln.^x = E Ih'" 



elfek 



fcGZ 

and the Banach space 

'^n,a — {h '■ Rk X Tg. — >■ C; real-analytic, ||/i||n.cr < oo}. 

To obtain a solution of equation ([7^ in the domain i?^, we need to solve equations of the form 
C^h = g, where C^ is the operator defined in (j64p . To find a right-inverse of this operator in (Yn^o- let us 
consider ui the upper vertex of R^ and uq the left endpoint of R^ ■ Then, we define the operator Q,, as 

^e(M(«,r)=^^aM'''(«K'^ (82) 

fcez 

where its Fourier coefficients are given by 



g,{hf^{v) = / e*- ("'-'") /jW (it;) du; if fc < 

Jvo 
^,(/l)W(t;) =- / 'e*7('"-^');^[fe](y;)dzi; iffc>0. 



The following lemma, which is proved in jGOSlOj (see Lemma 8.3 of this paper), gives some properties 
of this operator. 

Lemma 4.8. The operator Q^ in (|82p satisfies the following properties. 
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1. If h e Xn^a, then Ge{h) G Xn,a and 



Ge{h) <K\\h\l,, 



Moreover, if (h) = 0, 



g,ih) <Ke\\h\\n,a- 

2. If h G Pi^n,(T '^ith n > 1; then Ge{h) G ^n-i,(T o,nd 



Qe{h) 



< 



K 



n-\,cr y i — Oi 

In next proposition we obtain a solution of equation (j79p using a fixed point argument. 

Proposition 4.9. Let us consider the constant ki > defined in Theorem \4-.l\ and let us consider any 
K2 > Ki- There exists Eq > such that for any e G (0,eo) and a G (0,1) satisfying (j47p and that 
^■y_^ £^~^ is small enough, there exists a function C E Xiu defined in R^^ x T^ such that is a fixed point 
of the operator 



J{h) = QeJ{h) 



(83) 



where Q^ and J are the operators defined in (j82p and (|80p respectively. Furthermore, v +C{v,t) G i?„^ 
for (u, r) e i?K2 ^ ^o- and there exists a constant 62 > such that 



IICIIi,, < bMe 



^+1 Va + £ 
(I-a)i/2 



a - 



Proof. It is straightforward to see that J^ is weU defined from Xi^ to itself. We are going to prove that 
there exists a constant 62 > such that J" sends i3(&2|/^|g''^"'" fi_ ti/2 ) C A'l^g. to itself and is contractive 
there. 

Let us first consider J'(O). From the definition of J' in (|551). the definition of J' in ([SH]) . we have that 

jmv,T) = ee:^(0)(z;,r) = g, i-^— {d.Q%v,T)+dM''iv,T)) 

Using that d^Q" is a fixed point of the operator T in (1751) and that d^Q'^ is a fixed point of an analogous 
operator J" , we can split J'{0){v,t) = Bi{v,t) +B2{v,t) with 



B2{v,T)=g, 



\ 8 
/ cosh^ V 



J-"(9,Q")(z;,T)-J-"(0)(i;,T)+r'(5„Q^)(«,T)-r'(0)(i;,T) 



To bound Bi, it is enough to recall that, in the proof of Proposition 14.71 we have seen that ||^(0)||cr < 
K\iJ,\e^'^^. Then, taking into account ([751) and ([77|) . one can see that 



cosh V /— u 



(j-"(0)(«,r)+.F-^(0)(z;,r) 



_^ y/a + e 



1,^ (l-a)V2 



Moreover, taking into account that (^T ' (0)) = and applying Lemma IT51 we have that 



\Bi\\i,, < Kl^ile"-' 



(1 - a)i/2 



29 



For the second term, let us first point out tiiat 

T*id„Q*){v,T)-T*{0){v,T) = -dyG, 
Using Proposition H??] and ([75)) . one has that 



/cosh V 



[dyQ* 



cosh V 



{d,Q*{v,T)y 



< 



K\fi\^e^^+^ 



Therefore, analogously to Lemma 



one can easily see that 



d.gJ^iO.Q*iv,r)f 



< 



Using inequalities ([77)1 and ([75]) . 



cosh V (^=* 



(^*(a„Q*)(«,T)-^*(0)(i;,T) 



(v-p^)\v--p_Y- 



1 — a 



Then, using Lemma l4!8l one has that 



\B^h,a<K\p\\^^^^ 



(V^ + e) 



(l-a)3/2' 
Therefore, since g*?-! ^_ ' ^ 1, there exists a constant 62 > such that 



|J(o)||,^<flMk 



»?+i. 



(l-a)V2- 



To bound the Lipschitz constant, it is enough to apply the mean value theorem, use the bounds of duQ^ 
and 92 Q"'"* given in Proposition 14.71 and Lemma IT51 to see that 



77-1 



Lip < K\p\e 



1-a 



Then, using that e''^^^^^ < 1, the operator J is contractive from S(fo2|/i|e''^^7T^W2 ) C Xi^a to 
itself and it has a unique fixed point C. Finally, to obtain a bound for dyC it is enough to apply Cauchy 
estimates reducing slightly the domain and renaming &2 if necessary. D 



Proof of Theorem \4.3\ Once we have proved Proposition l4.9l it only remains to obtain the inverse change 
given by the function V, which is straightforward using a fixed point argument. D 

4.2 Proof of Theorem [2Jl)] 



The first statement of Theorem 12.101 is a direct consequence of Corollary 14.51 taking a = 1 — Ce'' with 
r <E (0,2). Note that the condition g')~i ^+^" becomes 77 > r + 1. The proof of the second and third 
statements, which correspond to r > 2, are considerably simpler, since we do not need to prove any 
exponential smallness. The first observation is that Imp_ ^ e''/^ < g. Even though in this case it is not 
necessary, we keep the analyticity properties of the parametcrizations of the invariant manifolds and we 
work in the domains D" x Tg. and D* x T^ (see ((55|) ). 

Theorem 4.10. Let us fix ki > 0. Then, there exists Eq > such that for e e (O,eo)i a = 1 — Ce^ with 
C > and r > 2, p € B{po), if rj — 3r/2 > 0, the Hamilton- J acohi equation (|50p has a unique (modulo 
an additive constant) real- analytic solution in 2?"^ x Tg, satisfying the asymptotic condition (j53p . 
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Moreover, there exists a real constant 64 > independent of e and /i, such that for {u, r) G 2?^^^ x T^-, 

\d.^T"{u,T)^duTo{u)\ < bilfile'^-^-/^ 
Furthermore, for (u,t) € I?"^ x T^, the generating function T" satisfies that 

\duT"{u,T) - a„To(«) - fie^M^iu,T)\ < 641/ipe"'^''' (84) 

where A4" is the function defined in (|57p . 
Proof. The proof follows the same lines of Theorem 14.11 We use the modified norm 

||/i|| = sup |e"2"/i(u)|+ sup \h{u)\. (85) 

uev^nlRcuK-u} ue'D"n{Rcu>-u} 

Then, using ([75]) . one can bound J-'(O) as 

Proceeding as before, it is straightforward to see that J^ is contractive from the ball B{fj,e^^~^^^^) to itself 
with Lipschitz constant satisfying 

Lip < |yU|£''-3'^/^ 

which gives the desired result. D 

The function T* satisfies the same properties in the symmetric domain P*^ . Note that D'^_^ HD^^ nM 
is an interval of size 0{e^ /2) centered at u = ln(l + y/2). Therefore, 1?^^^ n Df^^ n R does not contain a 
fundamental domain. Nevertheless, it suffices to deal with this domain to compute the distance between 
the manifolds in the section x = 37r/2. 

From Theorem 14.101 the formula of the distance follows. Thus, to finish the proof of Theorem 12. 10[ 
it is enough to use Corollarv l2.6[ for r > 2, to check that the Melnikov function A^(ln(l + \/2),t) has 
simple zeros so that the manifolds intersect transversally. 



5 The quasiperiodic case: Proof of Theorems 13.51 and 13.7 



As we did in the periodic case, the proof in next section includes at the same time the results in Theorem 
13.51 and in the first statement of Theorem 13.71 Then, Section 15.21 contains the proof of the second and 
third statement of Theorem 13.71 



5.1 Proof of Theorem 13.51 

We follow the same approach as in the periodic case. Therefore, we only point out the main differences 
with respect to it. We perform the symplectic change of variables (05]) to Hamiltonian p^ . In the new 
variables, it reads 

eKiu,w,9i,92,h,l2)=eK(xa{u),^^,ei,92,Ii,l2)- (86) 

V yom ) 

As in the periodic case, we look for the perturbed invariant manifolds as graphs of the gradient of 
generating functions T'^'^iu, 01,62), which are solutions of the Hamilton- Jacobi equation 

Kiu,duT,ei,e2,dg,T,dg,T)^0. (87) 

This equation reads 

e-'dg,T + e-'jdg,T + ^^^ {d^Tf ^ + ^le^^{u)F{B,,e2) = 0, 

o cosh u 
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where ^(w) is the function defined in (15^ . 

We impose the same asymptotic conditions (|53)) and ((54)) as in the periodic case. For fi = 0, the 
solution of equation ([55)1 is ([55]) . We study the existence of solutions of equation ([55)) in complex domains 
which have points close to the singularities u = p-,'p_. Nevertheless, in this case, we only need to stay 
at a distance of order 0{y/e) of the singularity instead of 0{e) as happened in the periodic case (see 
|DGJS97] ). To this end, we define the modified domains 

£>" = {u G C; |Imu| < -tan/Ji (Reu- Rcp_) +Im/9_ - K^/e} 

(89) 
Z?^ = {u e C; |Imu| < tan/3i (Row- Rep_) +Imp_ - Ky/e} . 

Note that the only difference with respect to the domains ([55]) is the change of e by -^/e. 

In the quasiperiodic case, the complex domain of the angular variables plays a crucial role. As it 
was done in [DGJS97I ISauOl] , we will prove the existence of these generating functions in very concrete 
domains in (6*1, 6*2). To this end let us define the complexified torus for any a = ((Ti,(T2) G R^ with 
ai,a2 > 0, 

Tl = [{6^92) e (C/Zf : \lm9,\ < a,] , 

with 

cTi =n- di^/e, (90) 

where r^ are the constants defined in pO[) and di > are any constants independent of e. 

The next theorem gives the existence of the invariant manifolds in the domains D* x T^ with * = u, s 
(see (|89p ). We state the results for the unstable invariant manifold. The stable one has analogous 
properties. 

Theorem 5.1. Let us fix any ki, di, ^2 > 0. Then, there exists Eq > such that for e € (0, Eq), a € (0, 1), 
H G B{ijlq), satisfying (|47|) and that g'?-! ''^^_^" is small enough the Hamilton- Jacobi equation (jSSp has a 
unique (modulo an additive constant) real-analytic solution in Z?" x T^, with a defined in (j90p . satisfying 
the asymptotic condition (|53p . 

Moreover, there exists a real constant 65 > independent of e and fi, such that for {u,6i,92) G 



Dl, X Ti, 



|a„T"(u,0i,02)-a„ro(w)|< ^■'''''^'' ' 



(Vi+V^)(l-a) 



\dlT-{u,e,,e2)~dlT,{u)\< ^^l'^'""" 



£ 



|a„T"(M, 01, 02) -9«ro(ti)-/^£''A^"(«, 01,^2)1 < ;r' „,,2 ■ (92) 



(Ve + V")(l-")' 
Furthermore, if we define the half Melnikov function 

,A„, /, /, ^ . /"" sinh(u + s) cosh(u + s) ^ /^ s ^ 7s\ , ,„ , 

7W"(M,6ii,e'2) = -4 / — ^ ^ i '- 2F(6ii + -,6'2 + -^)ds, (91) 

J-00 (cosh (w + s) — 2Q;sinh(-u + s))" 
i/ie generating function T" satisfies that, for (m, 0i, ^2) G Z?" x T^, 

(l-a 

The proof of this theorem is deferred to Section 15.1.11 The parameterization of the stable manifold 
has analogous properties. In particular, we can define 

. ^,/ r. ^ \ . f^°° sinh(u + s) cosh(w + s) ^/ s ^ 7s\ , ,„„^ 

7W"(u,6ii,6'2) =4 / — ^ ^ ^ '- jF[9i + -,02 + — )ds, (93) 

"'0 (cosh^(u + s) - 2asinh(M + s)) ^ ^ ^^ 

and then, for {u,ei,e2) G Df^^ x T^, 

1, |,,|2p2r;-2 

\duT^u,ei,92) ^ d,,Toiu) - f,e^M'{u,ei,e2)\ < ^ ^. (94) 

(1 — a)^ 
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We consider the function 

Aiu,ei,02)^T'{u,ei,02)~T^iu,e,,92). (95) 

This function is defined in i^^i x T^, where R^ is the romboidal domain defined by 

R^^D^n Di, (96) 

where Z?* are the domains defined in ([M|) . 

Subtracting equation ([55]) for botli T" and T", one can easily see that A G KciCe for 

£e = e-'dg, + e-^jdg, + fSS^ {duT'{u,0,,e2) + duT'^{u,eue2))] du. (97) 



Since Theorem 15.11 ensures that the perturbed invariant manifolds are well approximated by the unper- 
turbed separatrix in the domains I?" x T^ and _D^ x T^, we know that the operator C^ is close to the 
constant coefficients operator 

C,=e-^de,+e-^-lde,+du (98) 

in the domain R^^ x T^- 

Any function which is defined in {u € C;Reu ~ a,Imu G [— ^'Oi^o]} x TT^i for any a G M, and belongs 
to the kernel of C^ is defined in all the strip {|Imw| < ro} x T^ and has exponentially small bounds for 
real values of the variables. This fact is summarized in the next lemma, whose proof follows the same 
lines as the one of Lemma 4.1 of |Sau01| . 

Lemma 5.2. Let us consider a function C(w, ^i, ^2) analytic in (u, 9i, 62) G {u G C : Rew = Rep_, |Imw| < 
Imp_ — Ky/e\ X T^, where a = ((Ji,cr2) with ai = ri — diy/e, which is solution of L^C, = 0. Then, Q can 
he extended analytically to {|Imu| < Imp- — K\/e} x T^ and its mean value 

1 p2tt p2iT 



(0 = 7:7^1/ / au,ei,e2)de,de2 

does not depend on u. Moreover, for k' > k and d'.^ > di, we define 

M= max _^ \duCiu,0i,92)\ (99) 

{u,6i,62)£\—ln\ p - -\- k' ^/e .Im. p - —H,'^/s\ xT^; 

where a' ~ (g\, g'.^ with <7[ = r^ — d!^\fe. Then, provided e is small enough, for (u, ^i, 62) G M x T^, 



|9.C(^,ei,e2)|<4Me~^^"^(^/^'"''-»^'^^ 

where c is the periodic function defined in p6p . 

To apply this lemma to study the difference between the invariant manifolds, following [SauOl] (see 
also |GOS10| ) we look for a change of variables which conjugates C^ in (|97p with C^ in ([M)) . 

Theorem 5.3. Let us consider the constant ki > defined in Theorem \4-.l\ and let us fix any K3 > K2 > 
Ki. Then, there exists eg > such that for e G (0,eo), ot G (0,1) and fi G B{fio) satisfying (j47p and 
that g')-! ^^_^" is small enough, there exists a real-analytic function C defined in R^^ x T^ such that the 
change 

iu,ei,e2) = {v + civ,ei,02),ei,62) (loo) 

conjugates the operators C^ and C^ defined in ()97p and (j98p . Moreover, for (v, 01,62) G Rk2 x To-, 
V + C(u, 01, ^2) G Rki and there exists a constant fog > such that 



\C{v,9,,e2)\ <be\^^\s''--- ^^ ^ \lne 



1-a 



\d,C{v,e,,e2)\<be\ti\s''-'^^-^^\\ne\. 

I — a 
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Furthermore, {u, 61,62) = (u + C(w, 0i, ^2), ^1, ^2) *s invertible and its inverse is of the form (u,0i,^2) = 
(u + V(w, 01, (^2), ^ij ^2) where V is a function defined for {u,6i,62) G R^s x TT^ and satisfies 

\Viu,6,,62)\ <h\fi\e'^-^ '^^ "^ llnel 

and that u + V{u,6i,62) £ Rk2 for (u, 6*1, 6*2) G i?^., x T^- 

The proof of this theorem is deferred to Section 15.1.21 

Next step is to prove the vahdity of the Mchiikov function. To this end, we bound (9„A(u,6'i, ^2) — 
A^(u,6'i, 6*2), where A and Ai are the functions defined in (j95[) and ([55]) respectively. As in the periodic 
case, as a first step, we bound dy {A{v + C{v, 61,62), 61, 6*2))— 7V{(w, 61,62) where C is the function obtained 
in Theorem [531 



Theorem 5.4. There exists £q and 67 > such that for any e S (0,eo); « S (0,1) and /i £ B{^o), 
satisfying (|47p and that e^~^ i-a ^-^ s'mo.ll enough, the following bound is satisfied 



\OAAiv+Civ,6M,0ue2))-Miv,6M\<'-^^^^^^ 

for V e i?„3 n R and (61,62) G T^, 

Proof. First we define the Melnikov potential L, namely a function such that 9„L = 7M (see [DGOO| ).As 
A G Ker£e, where C^ is the operator in (^T)). by Theorem 15. 3[ the function $(i;,0i,02) = A(f + 
C{v,6i,62),6i,62) ~ L{v,t) G Ker£e, where £e is the operator defined in ([M]) . Therefore we can apply 
Lemma l5.2l 

To this end, we have to bound dy^{v,6i,62) in the domain i^^a x T^. Wc split $ as $(1;, 6'i, (^2) = 
'^i{v,6i,62) + <^2{v, 61,62) where 

^iC^-, 01, 02) = A(t- + C(t;, 61,62), 61,62) - A{v, 61,62) 
$2(«, 01, 02) = A{v, 01, 02) - L(u, 01, 02), 

where A is the function defined in ([M]) and L is the Melnikov potential. 

To bound 9«$i, one has to take into account that A = (T* — To) — (T" — Tq) and therefore, it is 
enough to consider the bounds obtained in Theorems 15.11 and 15.31 to obtain 

|9.$i(f,0i,02)|< 'f ,2 Il"-H 
(1 — a)^ 

for [v, 61,62) G i?K2 X T^. For the second term, it is enough to use bounds ((M)) and (|M|) to obtain 

1^.^2(^,01,02)1 < 'f „ llnel 
(1 — aj^ 

for (v, 01,02) Gi?«, xT2. 
Therefore, we have that 

\dMv.ei,62)\< 'f %, llnsi 

for (v,0i,02) G i?K, X T^ and then, it is enough to apply Lemma [5.21 to finish the proof of Theorem 
[Ql " D 

From this result and considering the inverse change {v, 0i, 02) = (u + V{u, 0i, 02), 0i, 02) obtained in 
Theorem l5.31 it is straightforward to obtain exponentially small bounds for 9„A(u, 0i, 02) — M.{u, 61,62) 
and its derivative. They are stated in the next corollary, whose proof is straightforward. 
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Corollary 5.5. Let us consider any K4 > K3. Then, there exists Sq and 67 > such that for any 
e G (0,£o) ind a G (0,1), satisfying (|47p and that g'?~i vg+v" jg small enough, the following bound is 
satisfied 

(1 — a)'^ 
forue i?K, nM anrf (6*1, 6*2) G T^. 



This corollary finishes the proof of Theorem 13.51 Note that when a is bounded away from 1, we just 
need the simpler condition e''~^(-\/£ + ^/a) small enough. 

5.1.1 The invariant manifolds: proof of Theorem 15.11 

Wc look for a solution of equation ([55]) satisfying the asymptotic condition (|55)) as a perturbation of Tq 
in (j55|) . As in the periodic case, we define 



Q(u,01,02)=T(«,(?i,(?2)-ro(u), 



which is solution of 



CeQ^F{duQ.u,ei,e2), 

where C^ is the operator defined in (|98)) and 



F{Ku,ei,e2) 



cosh^ u 



h^^^ie"^{u)F{el,e2), 



(101) 
(102) 

(103) 



where ^(u) is the function defined in ([5^ . 

Wc devote the rest of the section to obtain a solution of equation (jl02p which is defined in I?" x T^ 
and satisfies the asymptotic condition ([55)1 . Recall that D" has been defined in ([55|l and cr in ([90]) . 

We use analogous norms as the ones in the periodic case. For analytic functions h : Z?" x T^ — > C, 
we define the Fourier norm 

where || • || is the norm defined in ((70)) . We consider the Banach space 



^K.cr = {h : 13" X T;; ^- C; real-analytic, \\h\\cr < 00} . 

First we solve the equation C^h = g, where C^ is the differential operator defined in 
operator is invcrtiblc in £^,0- and its inverse can be defined as 



ge{h){u,e)= I h{u + s,ei+ 6-^,02 + ie-h)dt. 

J —00 



(104) 
This 

(105) 



Lemma 5.6. The operator Q^ in (|105p satisfies the following properties. 

1. Q^ is linear from f^cr to itself, commutes with du and satisfies C^ o Q^ = Id. 

2. If h e £k a, then 

\mh)L<K\\h\\,. 

Furthermore, one can bound each Fourier coefficient Ql: {h) with k y^ as 

Ke 



Q^'^Kh) 



< 



\k-i 



/il'^1 



3. If h E £^,(7, then duGe{h) G E^^a and 



\duge{h)\\^<K\\h\\ 
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Proof. The proof is analogous to the proof of Lemma 5.5 of [GOSlOj 



D 



Wc can obtain solutions of equation (jl02p using a fixed point argument. Theorem 15.11 is a straight- 
forward consequence of the following proposition. 

Proposition 5.7. Let us fix ki > 0. There exists eg > such that for any e G (0,eo) and a G (0, 1) 
satisfying condition (|47|) and that g'j-i i_a *'S small enough, there exists a function Q defined in 
D" X T^ such that duQ G ^ki-o- ** a fixed point of the operator 



T (h) = d,,g,T{h), 



(106) 



where Q^ and J- are the operators defined in (|105p and (|103p respectively. Furthermore, there exists a 
constant 65 > swc/i that, 

\\^uQL<b,\^,\e^, 
||a^QL<65lMk"-i 

Moreover, if we consider the half Melnikov function defined in (j91|l . 



|5„Q-//e''X"t<55|M|V 



2 2r,-lV£J_V^ 



1-a 



(107) 



Proof. Let us consider kq < ki. It is straightforward to see that J- is well defined from £ko,o- to itself. 
We are going to prove that there exists a constant 65 > such that F sends B{hz\pL\£'^) C ^ko.ct to itself 
and is contractive there. 

Let us first consider J" (0). From the definition of F in (jl06p . the definition of F in (|103p and using 
Lennira l5.6[ we have that 

J"(0)(u, (?i, 02) = a„ee-^(0)(H, 01, 02) = -/^e^Se (*'(w)F(0i , ^2)) • 

To bound it, wc bound first each Fourier coefficient and we take advantage of the fact that {F) = 0. 
By Lemma [5. 6[ formula PH)) and recalling that ^'(u) = /3(u), where /3 is the function defined in (PU)) . 
satisfies 11 ^'11 < iiT, we have that 



-^ (0) 



[k] 



< 



K\fi\e^ 
\k-uj\ ' 



-ri|fci|-r2|fc2| 



and therefore we have that 



J-"(0) <K\^\e^+' J2 



\k-uj\ 



-di|fci|Ve-rf2|fc2|\/e 



fcez\{o} 
To bound this sum, we split it in two depending j/c • wj > 1/2 or |fc • ajj < 1/2. For the first one, 

E -*- ^-dilfcily/F-dalfcalx/F < 2 Y^ g-rfi|fci|v^-d2|fc2|V? < :^_ 
Ifc • fi)l ~ -^ ~ P ' 



fcez^{o} 

|fc-tj|>l/2 



fcGZ2\{o} 



For the second one, we take into account that, for a fixed k2, there exists only one ki satisfying \k-uj\ < 1/2. 
Moreover, it satisfies 

-\ki\<—f\k2\ + \k-u;\<-^\k2\ + ^. 
Therefore, with this inequality and taking into account that for any k £l? \ {0}, 
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for certain C > 0, we can sec that 



^-^ \k-uj\ - zl^ I I - £ ■ 

|fc-w|<l/2 

Then, one can easily see that there exists a constant 65 > such that 

cr Z 

To bound the Lipschitz constant, let us consider /ii,/i2 G i3(65|/^|£'') G £ko,o-- Then, one can proceed as 
in the periodic case recalling that now 



u — i7r/2 



u — p- 



< 1 



z7r/2 — p_ 



u — p- 



<1 + K 



\foi 



and 



Therefore, 



u — p_ 



< 



K 



\/T~~^' 



j^ (/i2) - ^ (/ii; 



< 



vf(Tb)(' + $iii"^+'""'*-'"i 



<if|p|e.-i4±^|H,_k,||,. 
1 — a 



Then, using that by hypothesis e'' ^ i-a i^ small enough, 

Lip7"^/f|Mk''-^ "^^^ <l/2 
1 — a 

and therefore J^ is contractive from the ball B{b^\p,\£^) C En^.a into itself, and it has a unique fixed 
point h* . Moreover, since it has exponential decay as Reu — > —00, we can take 



Q(u,r) 



h*{v,T)dv. 



To obtain the bound for i9^Q, it is enough to apply Cauchy estimates to the nested domains D" C -D" 
(see, for instance, [GOSlOj ). and rename 65 if necessary. 

Finally, to prove ([71)) . it is enough to point out that F (0) = pe^^M^ and consider the obtained bound 
for the Lipschitz constant. D 

5.1.2 Straightening the operator ££: proof of Theorem 15.31 

As we have explained in the periodic case, we just need to look for a function C solution of the equation 



cosh 7/ 



1. 



u=v+c(v,ei,e2) 
Taking into account the definition of Tg a-nd Q in (j55p and poip , this equation can be written as 

CeC = J{C) (108) 

where 



J{h){v, 61,62) 



cosh u 



iduQ%u,6,,e2) + duQ''iu,ei,62)) 



(109) 



u—v-\-h{v,6i,62) 
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To look for a solution of this equation, we start by defining some norms and Banach spaces. Given n € N 
and an analytic function h : R^ x T;; — > C, we define the Fourier norm 



■■'^ ~ z^ II 

feez2 



\ki\cri + \k2\(J2 



where || • ||„ is the norm defined in ((8T|) . Thus, we introduce the following Banach spaces 

Xn,a = {h : i?K X T^ ^- C; real-analytic, ||/i||n,CT < oo}. 



To find a right-inverse of the operator C^ in (j98p in i?^ x T;;. we consider ui the upper vertex of R^, 
and uq the left endpoint of i?^ • Then, we define the operator Q^ as 

ge{h){v,e)^ Y,Ge{h)^''Hv)e"'-\ (110) 

fe6Z2 

where its Fourier coefficients are given by 
g,(/l)M(^;)= r hio\s)ds 

Jva 

g,{h)^''\v) = - f 'e'^"(^-")/iW(s)ds iffc-a;>0. 



V 



The following lemma, which can be proved analogously to Lemma 8.3 of (GOSlOj . gives some properties 
of this operator. 

Lemma 5.8. The operator Q^ in (|110p satisfies if h €z Xi.a, then Qe{h) G Xi^.a o-nd 

g,{h) < A' 4= II /ill 1... 

0,cr vl — a 

Next, Theorem 15.31 is a straightforward consequence of the following proposition. 

Proposition 5.9. Let us consider the constant ki > defined in Theorem \ 5. 1\ and let us consider any 
K2 > Ki. There exists £o > such that for any e e (0,eo) o.'i^d a e (0,1) satisfying condition (|47p and 
that £^~^iz^ *s small enough, there exists a function C G Xiu defined in R^^ x T^ such that is a fixed 
point of the operator ^ 

j{h)^g,j{h), (111) 

where Q^ and J are the operators defined in (jllOp and ()109p respectively. Furthermore, v + C{v,9i,92) G 
Rki for {v,8i,02) G Rk2 X T^ o.''^d there exists a constant 65 > such that 

||C|lo,.<&6|MK'-^ ^^f' |ln£| 

||a.C||o^,<66|Mk"'^ '^^f' |ln£|. 

Proof. It is straightforward to see that J is well defined from A'lo- to itself. We are going to prove 
that there exists a constant &g > such that J sends i3(6g|/i|£^'~ 2 ^^^^" | ln£|) C Xi^a to itself and is 
contractive there. 

Let us first consider J{Q). From the definition of J in (jllip . the definition of J in (jl09p . we have 
that 

JiO)iv, 9,, 92) = GeJm'^^, 01,62) = g J —^ idvQ'iv, 9^,92) + d^Q^iv, 9^,02)) 
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Then, it is enough to apply Lemma 15.81 and Proposition 15.71 to see that there exists a eonstant be > 
such that 

II iiu.cr 2 1 — a 

To bound the Lipschitz constant, it is enough to apply the mean value theorem, use the bound of 9„Q 
of Proposition 15.71 and Lemma 15.81 to see that 

Lip<A>|e''-i^i±^|lne|. 
1 — a 

Then, using that g'j-i i^^ I In^l ^ Ij the operator J^ is contractive from i3(bg|/i|£^'~ 2 ^^^^" | ln£|) C 
A'l^CT to itself and it has a unique fixed point C. Finally, to obtain a bound for d^C it is enough to apply 
Cauchy estimates reducing slightly the domain and renaming bg if necessary. D 

Proof of Theorem \5.S\ . Once we have proved Proposition l5.91 it only remains to obtain the inverse change 
given by the function V, which is straightforward using a fixed point argument. D 

5.2 Proof of Theorem [37fl 

The first statement of Theorem 13.71 is a direct consequence of Corollary 15.51 taking a = 1 — Ce^ with 
r G (0,2). Note that the condition g')^i i_q becomes, as in the periodic cas, 77 > r + 1. The proof of 
the second and third statements, which correspond to r > 2, are considerably simpler, since we do not 
need to prove any exponential smallness. It follows the same lines as the proof of Theorem 12. 101 

Theorem 5.10. Let us fix ki > 0. Then, there exists eo > such that for e G (0, Eq), a = I — Ce^ with 
C > and r > 2, fi £ B{fiQ), if rj — 3r/2 > 0, the Hamilton- J acobi equation (|87p has a unique (modulo 
an additive constant) real-analytic solution in D"^ x T^ satisfying the asymptotic condition (j53|) . 

Moreover, there exists a real constant tig > independent of e and fi, such that for {u,6i,92) G 

\duT^{u,e^,d2)-duTo[u)\ < feslMle''-^'^/". 
Furthermore, for (11,61,62) G P"^ x T^, the generating function T" satisfies that 

|a„r"K0i,02)-a„To(u)-Me''Ai"(w,0i,02)| <;^8l/i|V''-3'-, (112) 

where M^ is the function defined in (|9ip . 

Proof. It follows the same lines of the proof of Theorem 14.101 We use the modified norm (|55)) and we 
bound F{Q) as 

\\F{Q)\\=KW^ (k+ r 1 _ ) < ^-^\^\e^-'^/'. 

\ J-u\v- P-\\v- P-V J 2 

Finally, it is straightforward to see that J^ is contractive from the ball B{^e^^'^^^'^) to itself with Lipschitz 
constant equal to 

Lip < |/.|£''-3'-/^ 

which gives the desired result. D 

The function T** satisfies the same properties in the symmetric domain I?* ^ . From this theorem, the 
formula of the distance d{6i, O2) follows. 
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A Some remarks on the singular periodic case 

We devote this section to give some remarks and conjectures about the singular case, namely when the 
Melnikov function does not predict correctly the splitting of separatrices. We restrict this discussion to 
the case < a < ao < 1 for any fixed ao, where Theorem 12.81 holds. 

The main point in the proof of the exponentially small splitting of separatrices is to give a good 
approximation of the invariant manifolds not only in the real line but in a complex strip wich reaches a 
neighborhood of order 0{e) of the singularities. In this paper, this is the result given in Theorem l4.1l To 
prove this theorem, one has to impose the condition e^~^{e + ^/a) small enough. Looking at the first order 
of the perturbed invariant manifolds given by the half Melnikov functions (|57p and ([55]) one can see that 
this condition is necessary. In fact, if £^~^{e + ^/a) is of order 1, the half Melnikov functions have the same 
size as the separatrix when u — p- ^ e. A more careful analysis shows that duTo{u) and the remainder 
duT[u, t) — duToiu) also become of the same size when u^p^ ^ e. This implies that, when e^^^{e+\/a) is 
not small, at a distance 0{e) of the singularities the unperturbed separatrix is not a good approximation of 
the perturbed invariant manifolds and then in these cases Melnikov fails to predict correctly the splitting 
of separatrices. The correct approach when e^~^{e + y/a) ^ 1 is to look for the first order of the invariant 
manifolds at a distance 0{e) of the singularities in a different way. As was first pointed out in |Laz84| in 
the study of the Standard Map, these new first orders are solutions of a different equation usually called 
inner equation, which is independent of e (see |Gel97b[[Gel00llGS01[lBS08llMSS10[lGGllp . For classical 
Hamiltonian Systems, the inner equation is a new Hamilton- Jacobi equation, which has been studied in 
several models in jOSS03[ IBal06[ IGOSlOi IBFGSll] . Nevertheless, in all these works the inner equation 
was considered for points at a distance 0{e) of the singularity of the unperturbed separatrix since it was 
this singularity which was giving the exponentially small coefficient in the splitting. Nevertheless, for 
system ^ one has to proceed more carefully since this coefficient depends on a (see Proposition l2.1l and 
Theorem 12. Sp . Therefore, we will obtain different inner equations depending on the relation between a 
and e. Namely, in some cases (0 < a < e^) we will have to study the inner equation close to u ~ m jl 
and in others (e^ ^ a < ao < 1) close to u = p_. 

We start with the case < a <C e^, which corresponds to a wide analyticity strip. As we have 
explained Proposition 12. 11 the exponential coefficient in the Melinkov function is given by the imaginary 
part of the singularity of the unperturbed separatrix. Namely, the analyticity strip is so wide that the 
Melnikov function behaves as in the entire case a = 0. Then, it can be easily seen that the singular change 
that one has to perform is given by m = «7r/2 + ez and (^"''*(z,t) = eT"'^(i7r/2 + £z,r) (see |OSS03j ). 
Now the function ip is the solution of a new Hamilton- Jacobi equation. If we let e — >■ in this equation 
we obtain the inner equation 

9rVo - -^ {dzfof + ^(1 - A* sin t) = 0, (113) 

8 z^ 

which does not depend either on e or a. 

Certain solutions lPq'^ of this equation are the candidates to be the first order of the functions lys"'*, 
namely the first order of the parameterizations of the invariant manifolds close to the singularity. Then, 
the study of their difference would give the first order of the difference between the invariant manifolds. 
This equation was already studied in |OSS03| using Resurgence Theory and in jBal06| using classical 
functional analysis techniques. 

The case a ~ e^ is the transition case. In Proposition l2.1l we have seen that the exponential coefficient 
is given by tt/2 but the residuums of both u = p± make a contribution to the Melnikov function. Let us 
assume that 

a(e) =a,e^ + 0{e^). 

The natural change to inner variable is given by u = iT:/2 + ez and the rescaling in the generating function 
by (p'^'^{z,t) = eT"''*(i7r/2 -I- ez,r), obtaining , taking into account (|2ip . a new inner equation 

dr^o- — (d^ifof + — + n — -sinr = 0. 

8 z"^ 2[z-' + 2ia^) 
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Finally we deal with the case a ^ e^. We first consider the case a independent of e, and from it we 
will deduce the case a ^ e^ with v € (0, 1). If we take a fixed a, we have seen that one has to study 
the parameterizations of the invariant manifolds close to u = /?_ and the singularity u ~ p^ does not 
play any role in the size of the splitting. Note that now the limiting case is rj = I, that is, we deal with 
equation 

e-'drT + £^5^ (a^T)' ^ + ^ie^{u) sinr = 0. 

8 cosh u 

As a first step, we can expand the parameterization of the invariant manifold T"''^(u, r) as a power series 
of £. It can be easily seen that 

fe>0 

where Tq corresponds to the separatrix ([55]). One can see that the terms of the series satisfy that for 
fc>0, 

duTk{u,T) - -t: 

(u — p-)*- 

and therefore they all become of the same size at a distance e of the singularity. Nevertheless, in this 
case one has to be more careful, since if one considers the asymptotic size of the power series terms of 
the generating function T instead of 9„T, we have that 

TqOu) ~ 1 and TkOu^r) '^ ._ for fc > 1. 

(u — p^p '■ 

Namely, at a distance of order ©(e) of u = p_ all the terms with fc > 1 become of the same order but 
To is still bigger. Therefore, it is more convenient to deal with the function Q = T — Tq. Now, one 
can consider the change to inner variable u = p_ + ez and the rescaling in the generating function by 
(j}{z, r) = sQ{p- + ez, r), obtaining the inner equation 

^ , , o , , coshV- .^ ,x2 '52(a) . 

ar(p + Ozcp-i [Oz(p) —p snir = 0. 



where 62(0:) is the funtion defined in 

Proceeding analogously, one can deduce the inner equation for the case a — a^.e'^ with v G (0,2). 
Recall that for this range of a the limiting case was 77 = 1 — j//2 (see Remark 1 2. 9 1) . Namely, we deal with 
the equation following equation in Q(u, r) = T(u, r) — Tq{u), 

1 2 
e-'drQ + duQ + ^^V^ (duTf + pe'-^^iu) sinr = 0. 

o 

As in the previous case, we study the inner equation close to m = /?_ and therefore the change to inner 
variable is still u = /?_ + ez. Nevertheless, now the rescaling in the generating function is given by 
(J){z,t) — e^~''/^Q(p_ +ez,T). Then, proceeding as before and taking into account the definition of (^2(0) 
in (|28p. one can obtain the following inner equation 

d^^ + d,cb-'-^ (d^q^f + —J^-^lsiuT = 0. 
4 {l + i)y/aoz 

One could expect that studying all these inner equations, one could obtain the true first asymptotic 
order of the difference between the perturbed invariant manifolds. 
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